(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 7.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     78585,       1486]
NotebookOptionsPosition[     77572,       1448]
NotebookOutlinePosition[     77928,       1464]
CellTagsIndexPosition[     77885,       1461]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"*", " ", 
     RowBox[{"ODE", ":", " ", 
      RowBox[{"u", " ", "u_x"}]}]}], " ", "=", " ", 
    RowBox[{
     RowBox[{"u_xx", " ", "+", " ", 
      RowBox[{"s", "\[IndentingNewLine]", "Exact", " ", 
       RowBox[{"solution", ":", " ", "u"}]}]}], " ", "=", " ", 
     RowBox[{"10", "*", "x", "*", 
      RowBox[{"(", 
       RowBox[{"1", "-", "x"}], ")"}], "*", "sin", 
      RowBox[{"(", 
       RowBox[{"a", "*", "x"}], ")"}]}]}]}], "\[IndentingNewLine]", "**)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"amin", "=", "1"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"amax", "=", "10"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"u", "[", 
     RowBox[{"x_", ",", "a_"}], "]"}], ":=", 
    RowBox[{"10", "*", "x", "*", 
     RowBox[{"(", 
      RowBox[{"1", "-", "x"}], ")"}], "*", 
     RowBox[{"Sin", "[", 
      RowBox[{"a", "*", "x"}], "]"}]}]}], "\[IndentingNewLine]", 
   RowBox[{"ux", "=", 
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"x", ",", "a"}], "]"}], ",", "x"}], "]"}]}], 
   "\[IndentingNewLine]", 
   RowBox[{"uxx", "=", 
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"x", ",", "a"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"x", ",", "2"}], "}"}]}], "]"}]}], "\[IndentingNewLine]", 
   RowBox[{"Needs", "[", "\"\<PlotLegends`\>\"", "]"}], "\[IndentingNewLine]", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"Evaluate", "@", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"u", "[", 
         RowBox[{"x", ",", "a"}], "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"a", ",", 
          RowBox[{"{", 
           RowBox[{"1", ",", "3", ",", "5", ",", "7", ",", "9"}], "}"}]}], 
         "}"}]}], "]"}]}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "1"}], "}"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", "Thick"}], ",", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"x", ",", 
        RowBox[{"u", "[", "x", "]"}]}], "}"}]}]}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Export", "[", 
    RowBox[{"\"\<conv_diff_u.pdf\>\"", ",", "%"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5021991167832737`*^9, 3.502199165625288*^9}, {
   3.502199445286195*^9, 3.502199449171276*^9}, {3.50219980583748*^9, 
   3.502199871761393*^9}, 3.502199941959586*^9, {3.50220002712248*^9, 
   3.5022000282904387`*^9}, 3.5022000861529913`*^9, {3.502240809997251*^9, 
   3.502240811235413*^9}, {3.5022408799051*^9, 3.502240907782466*^9}, {
   3.502241322595303*^9, 3.502241326553067*^9}, {3.502246169014613*^9, 
   3.5022462300123243`*^9}, {3.502246552293344*^9, 3.50224658778115*^9}, {
   3.504000165207118*^9, 3.504000168624481*^9}, {3.504000224676523*^9, 
   3.504000253612252*^9}, {3.504000609202627*^9, 3.5040006192147408`*^9}, {
   3.504001766440124*^9, 3.504001788243786*^9}, {3.504001851228669*^9, 
   3.504001863656321*^9}, {3.504001897560314*^9, 3.504001903437768*^9}, {
   3.504002018033091*^9, 3.504002031573604*^9}, {3.504002110573649*^9, 
   3.504002111222156*^9}, {3.504002154418977*^9, 3.5040021909516487`*^9}, {
   3.5040022954364433`*^9, 3.504002299409359*^9}, {3.5040023338397017`*^9, 
   3.5040023377593927`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"10", " ", "a", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "x"}], ")"}], " ", "x", " ", 
   RowBox[{"Cos", "[", 
    RowBox[{"a", " ", "x"}], "]"}]}], "+", 
  RowBox[{"10", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "x"}], ")"}], " ", 
   RowBox[{"Sin", "[", 
    RowBox[{"a", " ", "x"}], "]"}]}], "-", 
  RowBox[{"10", " ", "x", " ", 
   RowBox[{"Sin", "[", 
    RowBox[{"a", " ", "x"}], "]"}]}]}]], "Output",
 CellChangeTimes->{3.504351095911498*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"-", "20"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"a", " ", "x", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}], "+", 
     RowBox[{"Sin", "[", 
      RowBox[{"a", " ", "x"}], "]"}]}], ")"}]}], "+", 
  RowBox[{"10", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "x"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", "a", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}], "-", 
     RowBox[{
      SuperscriptBox["a", "2"], " ", "x", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{3.504351095914826*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], Thickness[Large], LineBox[CompressedData["
1:eJwV13lcTG0UB/AUbZIWoR0VlWSpiN76oSjaFSpUEhVpDyEpFa3aKe2LKG3a
ZA2VLa1aKC0zd6ZBGTS30PZef83n+7nP3Jlzn/Occ+5KJ899x7m5uLgOzuPi
+ve59zirvf7Lfr11Z2Nru5//+m++vpTXoKQmljZrTifGlulpZN1Mr5Q0gOqn
g73Xo1r1CBnF4LOrzcCycVflfjKsZ/y+536GpDXu7y7u2R/O1lOMli2Yk7dF
rbbKmd+Zv/Wm9x67cXS1PYytsvfVPZiHr6EfGo6mHMWKmdmNorL8OOg8+SRB
0hlKPNHGRhcXo0FfqvZl+gmoaBleFf0mjg0KuuUceTes580Y35S8HOncjneV
8k7Bz9Um6+RjWfDTQnIPrPZA3p96e9PylXg41PH+pbUXePD+WUe9IrwsnPce
TvGGw8D7yJrFyvi8/qp+jKQfyuNaHx5xXY+krOX1Sof8sSR7pmBe3UbsXVz0
39P0M8i+/qCV0a+Bmu/NWmz5AMgLpHs0BGvD3d7+/tWj5yFU1z5Rxq2DVS1s
9RV5F7Bc2pD3eq4uYkvElC1XX0LDR5/voVk7sEs2P/+LSxDSTgjZvZPVx6Wm
NJHn1pcx1W2eerHYAJqb1WJsU4JR09+SLUIaIlvCPDxSMhR3FCTWpmw3x4X3
NxoOPAyFYJ7dbka6BQ6EDXErHArDs5uxo0o5lhAivYMepYfjB321bGK9Fc59
SDo3Kh+BZSpit1l2NrCK/lzzoD4CXyyzJUOSbKFusJoMPRqJBA1t/souOxCV
td6yeVFQX3arp9PjCCwTP500Wx0L+z29jt8Dj0LNWOGu1KtYiBfZvp7VdgI/
j/sI0+U6tosane6cdIJD3WzFM+s4FJwNLW287AwJlTh1m5R43E3Nb6x85oJg
/so1EZJJWHp4/yf6Dw+oGsZfu++VBKkXz/lqpz3REeb5pe9VEiJom/p9wr2g
yKNWrH4mGZIFRaeUM7zRNJOv/qEjBT63dbZlEr4QGk/Rko9OhYaIpbx63zlU
b/RPMaKn4kX7l3fFVwJg72U16b0tDcgqzDRVP4+yscV1Daw0KJ0sNpuKuIB9
rGv/ndyVjgDVobkxu0u4+fm8fvVsJu7uu0S33hiCHTK2+QP7syBYXd0clxKC
r3ZbFvCXZOEV35XJj9MhkBsNOFayLhvxz+b6E5qvwNNtsxJXfzZMtt/0SzkX
BjHnsjsF2rn4Ia2qkcMdARvbnLIfP/ORoioVY+URB3aftWCuVAFEzJrdeUri
EG7Pf9xKvwB3/EbOKI7GofKYp1RNUgGcWlhRdafiIeyhG3Zhy214hz3TV/NK
QGPIRzvewEJ4uTTpLI9NwqZ7orwy/EUYiOdPVJNNxRv1Rsf3G4pQGBzl1uCW
CseKc48u2RbBc3uM35uaVMTUDHoNFxXhVPf0Hal9aWDVl/TdNi3GS9du0e6Y
W8js2lOxMfEeHip3iiZLZmLhXPBhQ7kyfI6lBxXV5+Cm20F7td1lyDewcX49
mgPFD2qOoh5lWOb7+NNyyVzo3e051vekDNHc63dY+uSi7m5Kda16OSJ25ZrL
rcnDmNXoY+OWcihbbDOSScvH/js33vkuug/jDJ2tegWFULL8znoZVYWjojo8
701KIRVkuM0lpwpxN7gthyJKIVKSHSVYW4XvysMlqq9K8ZffSt2SVoXFKZcN
nXeWoa2+1ndAuxq7Jbd6VMmW4/yG4NnfjGqMCm8265SuQMviJUvUt9diyrbZ
ZbNaJc62bMVNTh1ayEXnfW7VYlmV6MC8ec/QMnNhqsO/HjU0/ebW/15g6lig
jrhWE16czJC2HW1A3sKFQj0XmkGu1l1ccvoVFLZmv91+qh3cf4pFpJXeIuFT
T2ymZBc2zznWWJU244Df7cDgul40h2q2NKm24tCQxqMt7v2Qj/v25+73NjzS
IXdNGQ7h2+TgQNeGDgjlPxa4tIOG/K+P3grHdkIrbR2zz52AsGwhmyA+QMb7
3qvPcUxEu1WfyFbshka/aoxhLAsT24sJ16AeLObZ+eti+Vf4bWqoWvqyFz2h
iRuN80cxWNoQ1b7yE84puBq63v+OqZMfbqmc78OGrvaCo5Y/sHYH3+yJR/34
9mHttpF1P2GYu9KlYtkA8uJ1cvqFf+GsrMOq4cODMDA+S69o/wW10f0Le0yH
kHxcbrbk8jjEF//gb/g7hI0mKd25qhx0PIm8umdmCFNnDSLH1ThIcFfibZ0b
wqnUbfP013Mg9taO59P8YSzwPbBlQINyWOMMW3gYRRVVDxbociAylfZLWmEY
PTs9KpebcyA0sqvf13gYFg8W3Db24WD+s1tlCunD0LJSjWqp4uDTfwZGg5nD
0D5gHjFZw0HZw9GhtJxhvAtN5cjVcWBXoysmVjiMO5djRdyecFBRMujHdX8Y
yj+Szb82cmCfrrDt86thdIsrnM3u5uBBQHFDyvgwZn3CenomOIidtDqyb2IY
fw9vvPj+NwfH/KfJRX+GEZtTGPj8LwfC3qZrwmaHodn8LSl/loPjruwIXwEa
Cn/l8B7kJSF2UMPcQp6GnHzvMHcJEqc0H30UMKZBQWCF24FNJOrbuu49MKVB
Zf31YGFNEhKnfwS5WNAQ7+o+0aBF4lmB0urG/TQYd0SPq20lIb4szueyIw3v
2lXVv4DEoz/OCyf9afjqqaUjZkpi4TNhXSKbhgVjZicVTpBwPKQikphHw0DK
vDW3XUhUT+rTd9ym4RhTeuMaNxIOGwIisopp0AssWqvoTqIyh95lV0NDOWdN
iYAPiUOhDzza39Gw6MeNgJBAEsV7nLKfTtDg823y75pEEppqK/Xwh4aK6Pk4
kkTiqfBQ37MpGnZ6KkzHJZNo6zyy7DkXHSK/pYbGb5DgHLGNfSlIx9GUawaF
6ST0fCwuvpajwy8Bag0F1Po0PZuO3XSErQob6qwhYRs4Q+7bQ8f+Z6uv99SS
oDk8Tuw0pkMqkr/44wPqforbWj9Y0HHBU0Ll40MSkmVau3vs6NjlvUSz8SkJ
p5dqmv0edJT6v9ewbCJBfpMSYabQ8exqo+nmDySSMkZepafSUWi29sLyLhIa
5lVBVul0mPNuev6bsvd9U3Z9Dh3ae4yTK3tIjJ273HLrHh0/n7daiPWRYCxg
Rlu+oPzpVKnvMInQ2vsG/I10WPG2q2+hkVBwC5p++ooOnYKQ+t+UjzYvd1d7
TwfH6bfEGYLE5wRjY75eOi5eO33RdoREl3yFwJMxOmzS77e/GCXh1x743PcH
HUnfjTe4jFH7f2VvgOo4HZV9YzcFvpOwZNJZN35T8fZdDDdik3h/b+lrHx4C
8Z6KTaU/SWp+oQep8BIQVXBw1f9F5YdI+ZYhfgIL2nlFuynv8d1TaCpMYDC6
zZccJ9G49WK4siSBayZ1PjITJJy/GekNShOIU3qlk0OZJ0NiIlmOQF1CEZ/i
JIkdXGXHeRQJvHyw5rbCbxJPmoYNBtQJjPX4/xT8S+WblSFPkgGBM+c8jgTM
kLhm/CBMxpCAcnX5+R7Kh/VV+Av2EJgbF+vWmCXBrbFwUY0ZAfa8pl0jlM3E
Wpf22hB4Ptsyq801gZWC21MdDxFA/B/uCMqceRXSX44QUI9x2tlD+dbPhBV/
nQgkjwV5e8ybAKvtgKrMaSoe5/7wq9wTePT61b18TwLDf8fs3lCOrddev86H
gPkdBxcBnglolktp6p0lwC9JrgmlzHcnqqYpgEA442jfU8qfsqa1zS8SuN5X
3zhJOfj6gK5jMIGKkRFd5/kTsL5qXs+6QiDC68/bG5SVg+p3eocTOMUfF/OG
cotHrmFIFBV//2DFmgUTyDkh/lYwlkAUhyZmTdnPPtQkMY4Ad9ds4SXKUmYn
LPOTCUxqS5xuoTy2q6dT7Sb1vLjt0scp1+saHahOI6C1ymR2Ke8EErXqenUz
qP3ui0vQpnxineqhpiwCVkHvbWwob1W69dksl4rPo9DiDGUhWSHHnnwCkeIN
AQmUK4W+O7PuEvh95pB9I+Xw+Q4jXvcIaFuslOynbDfd6vanlABjaT3PT8rr
ONtHgysI8J3/LDufbwJcoxUeglUEOFryLhKUO+mrfibUEHBesnZQkfLtvkRf
6ToC7qO3wzZRDuicP5H3iMDFXG07Pcom7/zPqT2lriuHHzaiLP+S+beqnoB/
oFa0BeVfDw8G6r4ksKV0lHWActP913ONjQQefLD0OUQ5tWhriNlrAi7LFyrb
U3bPLZrf85bAitKnAg6Ut6dJX3V4T0CtV2zZv+viCdECrFYC+swS83/fZ0bM
RHl1UPmutbny3/3rgj2E/3wgULZh345/vx8dMBgX3EPtt/mdP4aUHbwtxAU/
EZBbPtinS3mT2/PkhH4CUzIZIxspLzi6abn0IIEXs9dk/sXfa5OXljdMwPPU
nvNLKN+zWCKrRhA4JxfGxUM5yCgsq4pJoHeIVsqmnm+HvevhDV+ofGMYcz5R
VvQ3lrz3jcC4f5JvA+WzUerda74TeFfsqnaP8tsc0cS8HwQcSvzk/u2nzAOO
ufw4Aa8S5z3/9tuzpUfoFklg7zfi3r98WDKVEZ7wl8CSw4Yyy//li2iwvvAM
lU8mIWtIKr/q1jhzRc4RCObMurVRtrdWPR88n4Fxb8/EYMp3SmpOewkz8HZm
kWcLld9/X6aqjoowIK21ZW0GZdNPF0dcxBmY8uWTOUn5F6/+UYflDGzJskib
pc7Pf46t1marGLgR480lRjn2zH3Rt4oMaEQ9temgzudwdHLLrjUMrOCVpsdR
Dq87ZKSrxgA7szaan3KbGEtHbTMDn83OqNKo879K5d3vQm0GfrRpCSZR9kdp
tYIOAyOeXpIGlKXc/dZLb2cgsHFRY+YcVd8auRQE9zJgPPDVUI+qNzV99MEw
Ewbuz4SmDVP1SeBXU/o8cwb6eAWkrlAulYtZ+seKAeH1VzPrp0n8PispyDrC
QKmwW5XSFIm9sdNNxxwZiLRzOP6Qqn8Z+YNXBp0YeBhTsd2U8s72gpluFwbE
0xX8T/8hEa266WejNwM8jQeab1D1dOXnvb15YQzQOpmjiVQ9jnyhU7b/GgMH
+cSOzVH1+lehWjhfFANqUTuLXCm/8BHWPBXHwAmJZUWaVL135u+4vukWA2ZS
xs73qf5wR8PW8Hk5A50L57kZfCEhIrlXzreSgblAy2fRLBIBs9tIxRoGFhbr
DnVQ/WfvG5m8a48Y2GVoedmGSeKb/dCseRMDir93qO+hk1CPdKkZ6KPi1XBl
PesnkeJpExM3wEA1/d0DguqHc9Z7nHcOM8B652vKR7lNfq3YbSYD+UPTz3d9
pPpp9ffTp38y0HboR3Mh1V+rhvyVpnmZGE89L8bVQmLbltAkyU1MDOi3Zs8+
JtHtJLk5Q5OJM0mCrd8fkfCNLe1ZsYWJAr/HXJ+pfn+P0Sul/B8T0xX6pyup
eUAuaV3O5t1MnOpZ6WRURfWTn90lVnZMiD618RcrJvGuSOVVTAgTb63OLFWi
5hHX7qeuomFMJMz5S3ZR88p8bmvB5KvU/0lbsPQKNc/o2V4yTY9mgn3l1t/e
eBLlfJ2dRSmUz0i4uUdT84PzxaFXRUwExS+8sTGY6mdyrX/mdTIRb4LlGtR8
JSzpU/2oi4mDSi7DQdT89Vxcwtu/l4l73xm6b4+TWC1wiMX6zMSLXKeyg04k
2OPM7lYWE4tXVLvtP0Qi5M1sZcYstb6OYPaZUPvpp+6xVXUEJ4brHcXUSdh5
tquMq42gh3b5o4gaiUUn/Rj31o8gsyjuo5AqCR+Hh4dXaI2g/JeCx4wSCZ29
u435to9gJuPxQL0siRZ5e+Wu/SPg5RaPahOi5qO3MTSv4BF05LsFiXzhoH9h
l+v0lRE8P/Arksnk4IWJDPtq+Aimp7OdHxIcxLUUTWVGjSD5mH7woSEOVDte
LXmfPAKJ6hMHr/RQ8/XHeYYqRSOY6CvWcqPm7aYR/+LhjhE8MbpeEZXFwc35
9n77FFlwWVWRLG3KweyxZVrtq1n4ssEswWsvB84v20hzFRa8GYZHXxpysCFE
/6ypOgvH2wz8nXZy8HpO+byRNgtdAyIe4Vs4mPw7HqRnwkKVTajr8ZUcHPgZ
EaXqx0JIZWB54Dj1vjJQncvdwAL/5AYdgfhxxAoFBtY2UeuluT5GxoxDUMfA
xv0NCzvTL3QKRI6D+2bHou4WFghnzc9zIeP4Zck+d/cjCzezPm1r9htHW6Oy
uQWbhXXXi3sZB8cRXXJrKkPqC5JPGHlJyYxjQWCw9VbvL3DeW0OcTf8Frrv2
ZTSxrxjkv+IqHv8TPmllMWIvviLlg2JvtP4PPFl7jV/U4xu2butKTc4Yg+he
H7UQ4VFAjKc2JPEr+te1dQY8HkWrt22a0CgT5ZlQJOzHsMlF64XcWjqOlG0W
lZ0dw5bj9meKrQZwOSwnuL/oO5ZEb36oy/yAEpbsLWETNuRiXqTdbWnCtYww
xwEzNm7GXwx4daYJzvvGlEot2dAuJra8lWuCzOPH5aYH2ajrtbW45NGIyOuH
mqKd2LCZ2N2bLdQA182pPwUD2OBrLo5s0X+OVaFLjPhus+F01+txjE4tZrUv
Luq5w0a0p4HebGgNPo7RO24Xs1FYUGii3VKNuIOVR3ZXsNErAwEBxyrMqe7z
DXvMxo3wRR9lgivQ3349k6eTjS+m/ylL+hTiQfjksc4uNm4FFh1XnS1Ako6D
Sl4vG5cMZ2JbIvJhXKBetXOAja9XQ9Nas3Lw8Nz7N8Ff2FAMLE5yjbmBlHVa
1y1H2Zi1ltCuy0yCNy3deiWbjaVjH6+tL4uHson7YD2HDZ6FSTpSt65i/rwP
+XGTbOhLuJSxaoIxVK1z0vEvG+LHW+/Ozz2Lxyfz1m+YYWNuFeNqeoELbsgv
JOfmqPhef9KMFbqq9z9Axc6q
     "]]}, 
   {Hue[0.9060679774997897, 0.6, 0.6], Thickness[Large], 
    LineBox[CompressedData["
1:eJwV1nk8VG8bBnDmtKhEluoNpRSlJEkqyXUspUQREhVJi/pFSMqSpSxZIiXJ
nqWyJbKELMkS2VvsM2eGhJhpEanwnv6az/czM+dc57nv5z7PqpMXD51m8PHx
nebn4/v3qXd6sLViyFSjPNdUrFJGcNcsbQkH1jIVVI4cectNr9XYkhAV+3yZ
DpqSlktmPPms0S+1xueK3AEISr3U9Yn5rbG/sT03bpkJ6nSLd5emzMWakOWp
M9LmuL1xxy3VK6L4q2dz31rOEh63960QuieBYd/3VdaR1pDZ1xJbvVgGZqcm
Su8sO4VdTpEN8nFrUaUtUfg69gwslqee8ac2Qmn1rmdj0ucg5qoi3H5RGbGM
E2myyf/BokuvboeBKgQ415MOy9ljQO/+I5Obaiim2hpfmzjASquWvWgV4GB4
Su9YpCOaV8Z+yFighd5NAdq3ljmjUMe+RDxUFxEJ/6uQPXoZQ1pKxtr39KAn
nK5eFuuC6NSzVvPeGqCA27CVJ+2KyNZln5ZYH8IFS8vcAGs3WElqBRdrmkCm
iae4MtkdiutsOu/sOIzQLNF1RnKeECw5iHp7C+xenpIydNYLCE0pWp55DJ41
0YtemXjDu2POw1A+K6ioKtwyj/TB0fNSz45MWCNx8UH/oGW+aODKWGRcPAv3
xvtVh4t94Tx5mHl12BaH/SjG6qN+uPDj6cNSh/MQ/OnoVRLrDwNFM9GHCXa4
+j7i6oh0IMTchvdoaDvBOKS34EVFIASj/3fjiOIlKOrI/fS1DkKl1l53GWln
9D8vdFyeHAxBhQ1J5SIuMLrbdf6AXCgGmuWLxtTcoLB/dZpEbSj2ijDfnDR2
hwBx4fPA2TDM08oMdHb0gFXRdE65yW3YHbuz8FqZJxbL31Y8EhmOBaJ1lX8a
feAj8Hxt4LIIFAvzl8RtDsB63fCbuQ4RsFT+I7XndwDa/C4OdddGICH/uI1S
zU2sIRQyFF3uYf2YzeOas0GomUpRfN8WCU/9XzuZrbcg+CNyq3TIAyzf/B91
IfEO8jdfjtzb9wCbupVORKy4C0sH4wlHtWhMNnKaTOLvIntUuKhqMBoB/pqO
8UkRODR4U/387licLjx1c09ZJKJ63bTzp+MRW73F8IpmDDSlzFOYpglIve5t
97UjBsMW22YLZCVA9cCE9m/HWKwYcbXJ2pgIy79bBMqexOHiOVVZvp5EMB9O
p7uLJ0L0VPaT1O1JeDLH431vURKOmD/M/votBarpmur/3X8CXrfJ/CSJVHBd
kDv75xP4WwqcNtZORUzRCWdd4zQ8t7koURCRisJrnN/lIukQst/l577tEfz+
963XMzID1dc7LeZcewz1lfcUGOVPoZwpMkdKIB0/fn302FOdizrF6hONSumI
+3BAzGg6FydyrpZ4mqeDf8n/CjO2PcetApYDOz0dovJWOYYZzzFYkdX9yCAD
URZaaocj8xD/YV/O5ruZUDx9unmrewEWzPgc012RDaMDr/+6eBUj6pyZpcKe
bGQ0St6erCjGmvcKJ0Tss3El6IhPNaMEGmntNt2l2fglbr1gR0AJitIi8wsV
n+Fixom44LsvMWo88nJ/0zOoqbkkVheVwfTJ/beXFuZi4dEV3uy4V5A14g6+
Ds7Dm68TLmUBNZDw0lU7+zAP9itC4svyarAoKzF4fmEelnbx4j+xa/BbwFjR
iJOHjm793sfqtWipKLzE3J4Pv5e3t/75UQs3JZ/pX5/y0aK01LHhXB2ahMXF
FclC6Gh+VA53bMCVph2IGiuC8O5yx0yXVizNE2Hy85djMKieLN7agQKOdkOz
eiU2P1rmqrCVQuX5OEnzkSpYvcrZXx3Tj59yu4Sz7Goh4ff96fnJQTAmMxZJ
ytbj07qiu1+HRqE6c6LA+GkDPgtrdSYc+oYGX5WmmvXNyPgYt3NH+g9I3/4y
mcZtgWSpsvfvzeP4MsFiflBqg6HksqO7bvxCynBJvVDoO2xtlbzz5+tvCC1/
zOvvf4+AwSJxvq1TCDmXfyZxzUdAIGHk+qYZjJMZ/bZe7Ri1b8D0XT7SWbkq
b8nrDsw/FPlr71l+8ozKysTwNx04fDjiYOp5fvKIqkewYFMHTjYkls7Y8ZM7
1VROMjo74Nr6rvupMz/J0EoR5nI7EGSvvWDch58MM/I9Xy3RiXkS7SFrY/nJ
dAedlc5OnXgzk/NZpJmfZD2tCm5d1YXfDf3cbkUG2TI6a0R5bReMkn5Opm9m
kBUKu/UjFLpwQPwX/1UVBpmYXi1otq0LxQ+PzxFUY5AnH9Xc6tHvQnXE/WbJ
3QxyIO5N2OcrXTga+ain04JBjoY03J1q7ILybsWKt74M8s/59zHybt1YR+xZ
d7aNQZqoHg9X9OrGdOjstMD3DDKTf8B/i283ghu89dI/MkjLqAnHXaHd8IgT
5H7qYpCvqiX2GiV144rbrsw9fQwyYKX12NX6bhQ+MMrM+sEgxT6O6r+R6MEG
11QBYTGCvPDQRatRugcfjmmxxsQJsurCzLa2NT0Y2ny8vmMJQboQIqt7FHuQ
ZBDXEy1BkJ1KKpM8rR44LDL8PFeGIBOC3VKX/tcDxpnJPX5KBLlBc+70mZIe
PCy3b5qvT5DP0jPWdlT0gPIxWPXRgCBVxA2N9tX0QCglVzbhIEGqf45K2dDa
g1nSVzXXGxPkgVvr938d6IGd9qeiNRYE6dRpEOUq2oshqSSz9LMEWeQYoRxy
rhe1mu6ZP70JUqN7+7Ep+16UrtZ1drlOkJU6vX72zr24q/uBPXaDIOuXynUa
evXi6mVO5rA/QfaUFnotvteL7LO+k0UhBMk3v7shvqIXI+Elaz7fJ0jdpFVn
c5YykWc6HeuaSeeNEpKrkWJCsP3xkugsglwZ+qe/axUTpqaf3F48Jcjfrh9O
zlJgokhTncF7RpBZRoGWZiQTs2cZztHOJ0gx/m+mU7ZM5Cssyvcvpe830Ssu
as/EsG1ImV8ZQY6M1L+Tu8REJkcv+3o5QdZ0pBgZXmPCZtWb3S6vCNL12RGD
5HAmOjTtNHWqCZJlVaGzr5iJoe7pfe4NBNlgmkVYljPhX6Wcc6SRIF/sj650
qmJiTlXp3y1NBHlnmzMZ28TErIMKZH8zQe4WXqfO4zAhLVZ3acM7gswoC9ty
bwELRck4b9JJkFF5Ht/TFrGwpsGDWNRFkL7p53LKFrPg3x/iWUfbMlJ706A0
C/sm7yxT7SFIEfuJ9TtVWIjs3cv6ySTIK8utZNjHWLhoSQYL9xPkSUOlwz3W
LFRKSyk/on3gBn9Q+xkW3pl3Nqp9Iki5oeRvjQ70/Vb49x0fIMj3eYMVxb4s
JK1J8A0eJMiKz0Vj+YEsjOq5NYsM0XklgtflhLJww+GYUCTt694bbz+OYuHZ
g+azUcMEqaTvZBWRyUK3Z4m/7whBSnlp3w3LYeGY5JKLv2nPzRWvDSpggbE1
X89+lCCZSwsVfSpYiCk9032QS5AhfX+m7N6zUKyy9cRvHkFeXdK4+VwnC04b
DlJWXwnSZl/86VNMFnZLEsavaatlk40WgyzMqR4W9v1GkINufrG6f1lQO57x
aeQ7QQ4f3mD7kI9CcZTvL80fBHniqxhDmqDgZnB6KoL2/tWftkoIUFCZ/adH
ZYx+3peNzffnU/gwFFboQ1v1cMG5xQspCCwYudFAWyYwIG6RKAW7wdJvR3/S
9ZJx2BYqTuF1wJK7ibSFXh5pXbCUwrj8X9k+2pNc+dlzpSgkW21feXKcIFtM
Gi5MyVLonBQ6rTtBkHu4eXPc19F5rp4Nd6f9MiAu8dd6Cmnhw0+zaD8utn//
YxOFSN89pQt+0etpYmbvqEzBW3fB022074xCgKdCwfhKW5g1bc+VIurDOyiY
iunIZdMeK5r8cFadgtx/vV1ttM8bcy5+0qAgNdfMe4w2NVI/z0aTQrtB2hLx
SYI09X+eTGlTMJjTGLeZdr107C7LPRRqzEoWG9Ami3zbu/dSGKiz8zxDu+CQ
naP5fjqvDrP9Gm2FEdMF7QYUbBfNlblL+4WJfnmQIYXsRJOBx7S1S7UuaRhT
2Fo5R6OYdpPsjrXfTSlMPlkqX0/bInRTd+oRCh5C1x+30x4Ylw0zP0rh2WbN
Ig5tJysp7YWWFETaDS2/0J6uFZ2oOEHhzt6syG+0A5XmZTjbUHhqanzyJ23x
BzOW685QMDcmX43TTuQfF+2xpbBIyun5v+8Vzo/UhP1HoU348/bv//K2cdy0
7SnIiMcbjdDW2dmpOOFAYW5E5FQf7ebkZk76Jfr/3i3qnf/yCtZEWrpQ9HlG
X7jhX17nl3qirhRsfs+//PJf3t7c6Wp3ul9mC9mn057anZbr6klB3MXs9z3a
N58mnNnoQyGsY2iZ17+8SyMl2Dfoemu8aTj9L69XSFOEPwV7jIvq/cs7eP36
3kAKnIkrIxv+rWfRxeHsUAqX4wPNB+l6Nq06E28TTtfTes1oJW2LoGOHlkbQ
zx+wTjyGttPRfcWeDygsPnFTSof29GvYK8dSaBT/MSVOO1BBVWYgnsJjxU9O
fXT/JU7JBBmkULhidGOpK22FM8s0+B9TWH5TS5ekXdgk/D0vjYK6V/b82bSb
Ev6aS2XT/VFTph9A97uFwI+FLTl0vR2vvtWhPeAw9OpGHt2/YULMGXp/TGl+
lP9SREE0ala5Pe2b6Q298S8pLBDwviVDW1zsdfihcgrnJoS+vqP3n0J/9mRR
FQV+qRTHTbQt/APrbjZTeLvaybaO3u8DXK9r6m0U3oXpd56j7Wjmsvnrewq7
o14x59K+uc7mgVkX3T9jbx6r0/OioF79nNwnCod23VQJpOeNmPBXgarfFLaz
zfWffiHIttKggH1TFOpy+yQ3075zQXZO8wwFYnBPaw4970TrLYiuWWw4mfwY
zaTnoahf9RRPiI13my5+C/xMkIv+RH+XXM1G9oZIVRaH3v9pWx0fyrJxxF/o
Ikn79pEWntw6NnxtfoYmsAlSuGD26OaNbLwu8vEzo+h54ujwWXcb7Usthk97
CVLw8+6eS/vZEHlv5z7RTr9/7lEWkwZseDlJ7thGO0THvdPTkA378M2U80eC
XJD07GOgKRup91Z/GXpPkPOPS7YmWLExHDP/V3ErQQq8+1r99hIbEq9uCU3W
E+Ss8pjs1bFs6Kbt3ztQQpBd6jp7WfFsdLmn5jcWE2R28QgV/ZCNe+qL83OL
6PUv2CUq+piNr3vl/7tcSJA5WSxnvlw2RGf6i3tz6fdV7Gq13lo2Lva6rVRL
o/eXa0ZV5A82ar27masjCDJ0wvj4oXE2okpyhcru0PP98t+fCyfZaDZW7zQN
//f8Bmv9ptn4uCQenqEEedqWF3hpHgf6GzQscm/S62225aChNAeHpyoXPvEg
yP9USjrn7edAkWtb8fYkPc9bPmS+MOBA2/xqxyprglxs99XrrCEH4r1VGZet
CLI8VVau2pQDi50efkuO0fVfetvJ+wQH29YYCe8yJciSyVMLJi5zMN5/r3G5
Lr2e5UK7+hM5uBitG6kkT78P951MLBvnIGhr5nTOAINUUVilgUkODkXPJhP7
GWSZENVd/oeDa29G40M49Pn13fGlr/j6cJX8XmDFZJBjx81DX8/vQ5RgpAXv
A4PUcDL0eLOiD3Xc+tq8Kvr30RpH2vb0Qba5WaA7kUH+/CKxaCCyD10xnNpM
YwYZEfe5NvZBH7QXaHK5hgxyy8E8L+PYPrxSs5pSPMAgHXMNeBUP+8D00up4
spc+7171borJ7IPZknvhN3YxyE+zB0KMKvtwsM875oMcg/wgnTOvdLQP81NU
7b9O8JP5xrpEhE4/DAO80mUj+EmvvX4JeQP98Pciw92r+chVvXodyX6foKIX
JulQMwO1bb4Ry5QHsIMXEhSXOYVjK5on+d8N4PKGn+HW4X/ws/4Wx8HnM6Qr
vod8WzuJqFmWzofWDKKoapbF6ZZxiDHzkxhVg7BdG27V5TKG2dd8THY4DqEv
cYvrqv5v4EuzzOaIDmPNvCdH6vl4cIrOviVaSdv6f6HK3CGUbrgpIGL/BRde
J8TMrP4EET0nhetCI7D5u6/++TCFno0t71xfjoDdoPz8gGUnnsVjTb/lKEDd
V2jRbsMPtx9zra1HUfDim2e2fBtUzR5/6bUZxd3IZ+xY4TaUCAs/b7cdxXTX
l18Pu1tR483UfOs0CjsVgUOHnVvRffKaVa7/KNYska2zf9KCOWtLor2fjsJr
eXKr4rYmHM9WFVk+PQq5iE31IlZvUGDxS6iFj4szsw9bhCi+gdCcYsEbBBe+
IzvKFaZrUXFsl8CQABesg3fMWHG1kJmvM5MnxsWpRoWxhawaDJwyGtWX52L5
jQrLJNtq2EtcqPMw4WJjmbfelgeVqKneWLvJjAuf+u/LwuwqIe3Iq+KYc3Fv
naCamFYlWmqdKvZacaHbU1gc8eUVlC+7Foqf58K7KYvjSb7CRLN/aqYXF9dc
v7/l7y+Ht99Dn550LtQvPz4aMFiMhVfm6llncYFHE/cqEosRbWsnOpDNxbkm
u5065sV4vn9HMi+Pvl6SsM6Gt0UYEGmtYpRz8Wjj3WfuuS+wP55PQP4dF35r
N87SDCvAkkKr0Mt/uBD7eD5R5Xoukp9UH/49xcVqkfFv49q5UIreIO1F1z11
6+iNuXNyse/aRPbN2Txs2jLBFgvOgYd2WGuMMP29T0tbT9QzcJrLxCtX83D6
98bUu/5ZyBpcHiOkz8PQ6Z+alU8e4Wac3wnmAR7+GD06uUnlEU4dGpV9asQD
6avcYVuRCqmXL58ZmPGQfjxlrlZ3CoLCjtaEnORh9gGpo5ESybBVffBtvisP
a9eP75T0ToDOF77CLncelGr5hBvo5pJOtPVI9+RBUWIoWbovDh/nbZ+r58uD
RXv0nRz9WOzu/SgVGMZD8FfXH+FbHkDGV3zv3Ec8HCy+EdtseBfT2z0Wtj+h
86SKC0/E3EHnaF/bowwemPcS5d8OhuO22fPje3J4WPh+bZ17wG3MrD90ye8l
D/3FX8Kkc4PR0xoWT7zjoaSUDFsc44MX/hM27z7wEP7k0sXgs96I2Gkln9zB
w/aj5+b22Htif6pinhaTh+LUMNm/la4ovtpY5zPEg5W+29WPpQ6I3Lg1zGiE
R59H3V8nStrBkRNrsorHw+sExVt2dbZYp3+BVTFGX2+HwGTKgROYxf8+5fYE
D78U01cq6x8Blb/zPH2wx5zGzaLflxji5fnkTUpTPCS5BzXojwP3pRf8nJnh
YZ+veK3HiKXG/wEn4cQY
     "]]}, 
   {Hue[0.1421359549995791, 0.6, 0.6], Thickness[Large], 
    LineBox[CompressedData["
1:eJwV1nc8Ve8fAHBX+iaRjAZRdpTVsorPKbKSrSJkpShEmraMSigJSZGQGcrI
yIOMEFn3XruScc81LvcaDfE7v7/u6/16nXuez/N8xnPEnbzML7CzsbHNk9jY
/v9reIHWVYtbaU0L+lww8tmuyaEtfPWb0CFwNRMhiVT0aB1MTUp5L6QDEqa2
i6rlLK0xEamQmzLG0Jiq16Ip/x+cbKe+eyFkCWJq1d9n3guA1EPRzLXd1jAi
X5H8umkXrBg6JzrK2ANK2PuvR0YG6GG9DY4JjkCrExHYL6UIZ1yWP8YJucCf
CX7xmIKD0KAtXP4pxRWWXwUqZoyrgbKkZtHCbjeQlJecOlipBSnsDjnSry/D
sTIF5ainx4FzNDT9tIwnRD78epBfRBcqv3e3f7K8CkviJ62ezxnAVVMXQ9sE
b2i+ndDFL2AMw0qR2tFCvnDnTX5eTZoFxKfuqJU+dx0uJH7bnyFzGgx5c4/W
pNwAqVwf2tmes1A2++UwY/dt8NpgcnHDK3u4Ym//LtLxDsTbKLotfXcAiQ6G
othrP7iW9ta6OswJYgr4Zc1kAuEEM+2Uwn+ucEI0IwO/GARuXpdZ7T8vQmBT
8pY6y2CQ0B7fUdDrBodU5KOtE0Lg1AG9Sb45D0jbahLxQCgMRC6rll6NuAZ+
7YkNpyvDwPhyzob5Ol84Hf6dXfJcOLjXtpv//e8GcC96B1WlREC2Vd6v9IJb
cKs3/tb07vuQeO8ybf5KIFg8HC77UHsfNr60vHN9IQgUdWQWwxwfgIHoh9Cy
+mAYe1/uLfo6Cm7E5FyiOIaC2ZMBd2OZGKjVmxZ41BcO8iclc4SbY+A0LRv+
1EUA57orkxMXY+G2uC2PQFEknK9YLUaWj6BXM73JMPk+bJV7pHg24THYad2m
52RFQwjn+z33heKhxT1PZNb4CezVe3zv3dV4WAoUkWtbfALd4V74YHM8GG9a
jLZJjQepdfJ5ijeeglid+e2Hf55C078Mxd7uBHib+f6UeFcScLMSDu9++Ayk
05qfKa2+hNL91xP0fz6DpljjXVVfU8H+qsWyt0YyuD2h3FowS4PCGd6KBloy
lPccLxbY+wrMafeOup9IAer+F1Ib6emQNHxHu5R4T6M/i9cmOwuOiVhnjFil
go+PzBmW9Bug26iu5yxIBeF20VPkjDewa/q2c4FCGnzyQNu9s7LBy01Fmm0o
Dd5sN+j89SEX+F0KszPV0mHbyeTsL8KFcNb6VeHcfAboaDw43uxaAoxBS650
4UzwNzfT5y8pgQh7zgsW2pkgk9DxZoxUCu+dvYTL4jNhC0nLb1dqKWz21Az3
U82Cf1GFL799K4PG0H6b/wLewA9b40KHOxVwIJ/vPxHOXNilcXUqhQ9Bi2Kj
Q7tyLlg1pf0444XAofhWVaB1LtTVIsm4dgTRZd+u/sjNhb+/DtXWnagFWm3B
YNapPPBpWS+R818dvCQbFO9/kg9D7N1OB+LrYdNaiK3erkJ4/j7hkdZwIyS5
nbGX1y0ETVsDhWr+JpDqlXfg8ywE3+XFHh/9JtDKoToPfiyE0f1YxYOSJqjI
SSgtVywCd3EPFcfYZpixmK4+2UHYxNcxyqQFrLIT267xvAOxVlL52K8vIG02
S/sUVQJ8fXnLa9e6QThIT+PiqxJYSwsB+aRu2FKQFsVVXgL2FNNJ3+pu+MNp
oWg2WgLHAtt+2azvgc7a8msjaqUQ4bsSczOpB+4oh6z+Gi8FXvWAL8OfeqGD
V1BQESuHf7Ss7eWKVLjZoQ5JCxUgKU0LNjAZgu0lfCMkEgLLo+s9rZ1/Qtmo
9pevR+thBp4qecTgUO/+Yqf1dANMweP3pTgDFmU0eQs8mqHrmqL7yhIT2H/n
bdkp3QrGymdiHBSWQGXNoczi7RfIW9zKJbDpN3wJO9TRtPcrSPg1N9Gy/sLu
R1O/c2Y7YfCLj9OowirUOAblvGN1gljVraBUlVWwPyhgXfWrE2J2XdxzBlYh
haLx4Qt7F5isO7qxzHQVdog+uMHY1gWtoUsF+6+tAn+uLOsQdAHHJ3Mj0fJV
WN9wYRrFdoHjU4ycr7EGGU9/P/8c3wV8E1wWJcfXQOdi9MmuZ12w/cQYXmq4
BqFcpbmjr7uAW3C31RubNVg15XBb/6ELwKC4/vidNfg1nD5u+L0Lvhwu3Wj9
YQ2mlr+NkJW7YS+31q6l7WxY2cXYekOVbrjPZy83sZMNC6ZqZaEj3SDOKR3U
vZsN21b+wiNHtxv8psvY0/awYcdvnFvxt+2GHqviA5tV2bDkBYqQ1L1uML9c
EHLYig3Tn2+39PnWDWK37R6nx7BhAo4BapNj3VD01lBtRxwbNtIpL2JL7wbG
naiZ+0/ZMN+iqJ8nFrvh1iuGn2MKG5bqZeAjtKkHhAt/L33PYcOWZhpia1V6
gJMz4sfuBjYsg17VujmmB9qTNtY+XmTDPBNNJN2f9ICBb+KOW7/YMFWdn36N
RB05Sj2/avOXDWt9waXg/7oHNs9KCAmQSNi8qXUs/qEHpgTdJS25SZhW+ZJF
w88e2DjBMSQgQcL67+4fvq3eC+HaFkrBhiQsXbnxMFmrFzy1zD/+NCJhV4bP
xijr9IKd08IubRMSxqYaDJPGhFuEzi5akDAZescrS5de0KVuDtpjR8J8Ta+4
KsX2ghnzZkydJwnbLPqGMTbWC3+FxsIaYkgYdvru8E+8F1K90MngRyTMJ/Z8
2+hsL9QInTJTiyNhFNKON99/9YKr6q6tL56SsJcT9+yHuMmwrUT9knoKCVMs
cm/vOUQGx2dN4dXZJMwB163qVifDkr80+UguCYuTkMzp0iKD1kuj0+V5JGwp
fjDsqz4Z2icm6elvSVjNbaOjbbZkaBg4qW9UQsJOaSvm1YeRwV73n9JkDQkL
9Od6VnefDGXBBtpytSSsqHQiojaGDKDa5n+xjoQJyqY61SSRQSql4nzfJxI2
xL1FuDKfDJ9fJHZHfSb2qzu9oaKYDIIvBjkKW4j9Bn1eLC8jQ7UJz9mOVhKW
MR/cVVpLhtM5hfakduJ8KfP3invJsHXP9UG1LmK/vB03ivrJcFAzfZ1mNwnr
1M91KRwhg9jGL3qaPSTsUJUTVkAjw8bt5eIHyCRsJbVnOXuFDMWJaZtm+kjY
Q7dS1zQpClifnVB/842E0bGCe5v2UkDJUvJSwHcSprcjM/eGEgU0DdV/Gf8g
YezN8bNGGhQ4HmVzZGyUhN2W8r3x24QC5SqNR6jjRD7+Xk5ysaLAt+e8F+9N
EOv3OFd+taHA4w3qjEOTJGwuxGI18wIFDv37KhBMI2Gu3w5EmPtRAGdufDlC
J2ENZXuzPwZT4DpVI819ioRJxEi0ykZQ4MjB1j/zhIeP8m9ee0QBUzFaGXOa
hGkIcim7J1DgpI3H9iszJCxpimROfk6BDXZaHd8JWybPP83LosD478aj1bMk
7J03Xr4tnwImvqEsCQYJ4zX40R9STAH1WwmrYYRblztFz1ZTYEf7OO+RORIm
+/UzfKqjwJSaytZYwhFZtY6KzRTY58LpMUL4mFVRJkc3BWwlvPm95ol8yGc3
e1EpcIlH/nwR4ZV1afjAEAWwS8Ur04RtBhM36Y5SoEAiiS7FJGEf3sUqFE9S
4IbgbdmzhLc9iDQRmaGA37/Z0gjC1xyDvCOZFKgKT4grItypdvMJc5kCrpqb
a8iEFbd4ldr9I84vZUl1ifDDSVfqZ3YqcKxnrvGxiHzW2P8+yEmFppqMrXKE
9RNO70zloULu35JbRwhnehhrcglQ4Rm5VdaAMPsJ3fPXd1BhyfiRtDlhBxGt
kO+iVMgxL/Y8TbiGdfj1SUkq7CsZXPd/i7QpNJbJUuGeEhozJXwnXXpSnLiH
PN14+PQIU2+Lbow+SDzPuH1PjfAhs637fqlR4ZNKpaE04ThZnlPOWlT4+TL8
DA/huTUOrw5tYn3Pp2/niP0YU1ceqRtQQYyaYfKVcP7bhXcZxlTQnbuqkUN4
Y8R0L68lFQSrP3gGEna1G1u6Y02FBSnNaWPCDYeGdkzYU0F8tfadMGEJ7l4N
MxcqJKzbWP+DyEfwzzbbajcq7B4YFMogPFL5KXCPFxVMeNkrHQkfiatKi/Ol
QoCQbrow4SUsb+xSEBXk+M8ZBhL5r86Yen4zjArLZPdhRcKhnPIWEfepEHw7
0mKQqKfNnfl1r59QIUrUmraHcO/BmVvvkqgAKa3/fSbqMTlRQbnuBRXWLddt
cCEsc/7ty5E3VGiQbyuKJep3un7WaiafCsMimWd3E34no8SzUkwFdYGno7lE
/WvNFvoJVRP523cmpZToj9MBxWctO6nwo2tr5jOcyM+PeV5nMhX4E/XOcRIe
1TnQ7D1ABWGNRaYP0X+e3O8Px45RITMirukI0Z8Rz0v4W39RQfNo27/nYyTM
aG2hpe8fFRo7m2vwnySM3/lwyCR7H1RUsjseJPxybxmDg6cPtE4oG1QS86G0
ovwLiPeBoURGdsIIUQ+iv+4ay/SBpfKna5+GifkWonbEbl8fKDwBsekhEvZF
vyLnzuE+6NL3E90/SMLGqJURpQZ98G7UL/c+lZifSx+xfT7E++n6u69/JWED
1qvL6jf74K38cU6HDhKW9lGrUN+/D76Sp3p0ifkoH45EXSP6QCOKKsDZRsJ0
BOv+pCX3wYSFgIZ1E3HfHGgo2fapDzQbcB2VaqJeDomlPf7cB/WFpKf1lSTs
rIp/FHdHH+zlyxAwrCDyr3HIib2/D17uzn1hWEb0x/EM3tnZPvist32Gr4iE
xZqFuTcK90N5dFt9RTpRTxbfrTCxfrhYkST+O424r6yOHquS7gcfP9uSg6lE
fqwXthcr98OuJxdvJj8n8uPo0pii2w/7As4GKcaTsNyrOmK+Pv1wdOba36th
JCzFJ20T42Y/NMv4dDiHkrBo35Ult4B+CBdsGTILJmFXb5W0n7/XD8XUFwW7
/UmYarCU38mX/XDJXkgp9BpR/7EcFInWfoieuLTf2ZGEfXvbENUlPgD1Ic7P
nqkT82SGY/rAngFYSDpnT1IlYbXyJ4zi5QegoCfqpcsh4nxzG7nPqA6A8f2A
J6JKJMwpqyl6yGgAWL2uLtZSJGzixefYyZsDgL6JG//mIWEzD788+dc+QCR3
9u+3YeJ7po17wb5nACqE5os/DLBhHVynLGv7BuD9v9Z/D6ls2Nt77YJ3fw4A
u1+HjWwXG3Y1rOMp5+8B0Ps4p7OX+J5Z8OtMFJAahC2nZcvOEd87f917n8vd
GQT9j0qG/lfZMEsVu8eKQYPw7IPr6Y0ebFg+aSLiYNggVK1JZMe5sWH2Scve
mjGDkPBiZneyExtW1yisb5Y+CIocjn7+lmxYpJjjwq3WQajlfVwZS3yvCVBm
jD4LD0Fl2s2G2cU12Hdsw6pr1RD08uno8s6uQlFu3p6+2iGo9Q7ISRtbhUOC
pmYGTUNgxLZhRm5wFY5OJmXs6xqCIaFqcanPq2Acvffk3ATx/69Hn5e/WgWf
/lNJt/mHgeUdl2lnsQoV3vEHHroNQ+7yphsfi/6BXrr4xeLtIxC/QVvhseUK
3BQ9L/HD9huIeDv3Owf+Bvlpq03UU9+B+0JkcL3jEgjwznE2/PkOf3a7V/x+
zwQO9LxQMuUHjFTxJGqIzMLlQ1X9G0+OQuLiW8nFiXHIM3BKq1kahRpasMUH
z0FYnBLeMpHwEygXdjqnhn6GUgu9dfE6Y8BTc2VjCedbFKQfnloyMQZdZuxX
Kpa6kPiwYd/r8HFwG8sS3TEzgjRUw+KFDkwAxz3hjvNfxxHFSUjlxaEJsD4W
dOze+Di6FvOWKqY6ARoN2u0Zf8dR/nifsOzRCRD1iuav3DOBdsUrvFLRnYA3
f92eWgZNIPZ5SoGFzQRotTx5oi0/idpy5ZqjQydgZCCebdCfhmx3ff1N6pkA
I6o3ur9MR5uFfEqryBPgMnGrS4R7CtUJbPW+3jcBH572caaJTSGZjedotOEJ
aC5War1hMIUYrAnKV9oETMUGtLokT6HQltX3L1YnwOKNzdWj6tMo21fRU33v
JLTuFtue4z6DbLy65Fjyk7AvT2TjwYAZxOPuO56vNAkxr/iO5cfOIJ/zlbZi
hyfhuPjbS1dKZtARQ92TG7BJ2Hwz+7vEvxnUsdtelmw1CS0ujv1OUbNosTV6
9GrIJCjUcgWve85AQ5vIl1buTgKXkoWfZS4D1RuJMCIjJuHLTNnhhAoGetSR
+/dl1CSoawN5ispAe7ubBdufTsLqY+1UxsY5ZN9P0pPLnQTR9GdO3jpzqGny
et6P7klQU6rkTn0/h/JlP+73IE9CwTaakknlHHrixvHhF3US5Gj3PJZr55DD
VFwD7/AkyJLUlJU75tCf2cJhzclJqLe8nXKcNoeUlnDe5L/E+7slnXHheZTE
Ye9rLkWDy3NC3fvvzKNV5+2Hu2RoIHKRS6ExaB65fOpcNJGjwZcxlQrT8Hmk
HKp985QiDbRzJM+bPppHn9dk7+ir0WD/37Wy+cx5tPyHFaRlRAP/TYqnqr7O
o9Pz96P2+tKge9nJ0kGUiapNtY1ybtBgazIH/7I4E0kUrXDL3qZB3H0OcrgM
E816esVIBxLPq9q4RikxUfi05WOxezTQ8nvyVOUYE5VM7k7clkIDzRMcctzO
TCQwUprO3kADQXO8LT+NiWK4AwLKm2igS2q5EpvBRFxHdM5eaaFBk3zy3UvZ
TMSe1M1D6aDBR8eGzauFTMQ0Y9zK6afBlO7fgLIaJupslDUxZdBARY0+em6A
iYwWGHLrmTS42b3ePXOYiT5LlHNULtDgXdubM+PfmQgFnaiU/EMD7LeCs/4k
E71Vd5Je5sBhWzNd6jmLiR4WPP/7QhiHuL+uBT5cLMQ55EQxF8VhpPD1MQke
Fgrj2lu8QQyHpVcM5c+8LOR38YOrtzQO8YNOu5YFWchdnNylo4zD++9W0j92
sdCESUr+7wM4HG6PDdUWZyGnQOfIt4dx6NFMUUmRZCHrwfmjO47gYHjT5pWy
LAvpP92cPXUCh4ah6OvjyizU8IkcmqaPQ+9Fx4A/B1gIY6bYWZ3EofhpW9/6
wyykarJPAJni8Pci//01NRbRb/qBcedwOFDN1xyIsVC6Kq+1nj0OH9aknXWP
s9AuV8rBFQcc9H4nmHLosNC2Ty74BVccoiZqdjnrsdD6gBBLdW8c5i+cvZt/
ioVK5qd/o2s4pKXFV/KYsJCL69lU3Rs4hHKUh7iaEvGZKNEt/HDo5PBAq+Ys
5NuQHDsQgMO/ksnDWpYsJKX+32HHYBz8y7YdvW5FnJ/EcJBnOA4rjSHLLWdY
6FCivvRiJA7JSR/zJ8+y0Nimkla/Bzic1i2YX7FmIZ3FB1vvx+KQWdmyd7Mt
C3F8Vst7/QwHP1bMubbzRHxHM0z3puDQQZ+KzXQg4ivmXSp6iYPssUeWNx2J
+JInsJrXOIRnHYwlORPx8ZpP6GThwN7gv76SsFTYx6i2bBz+2L3idndhobse
8ZS+AhwcdQ92FF1goQOja37ni3AIuLnuvp4rC42eviw+8Q4HFt2XTCYc10Zp
ulKKgzWl9v25iyx0HDt+hVWOw4YoSal+wsySAr47lTjo3BiSMb5E5ENOqJzt
Iw53xRWrKwmbvwyzjUQ4vD5gP7bLjYXYBeZIm+tx2HWlNPUO4XeR597EN+CA
FQTNthN2Wmky2tmMwwXmpq4d7izE732A+aoFB2mrfL1zhOvHXyTKfsHhIXeZ
+VPCPjYbNQs7cDjlEcJsIizx1Xf0cBcOi7jbnjnC3drfI6t7cBAW6JjfcpmF
Qj+cVNCm4BCbyGsuR/iAQnl3Sx8O0TpP9NUJj76SuGU6iMPU+joyRjhuW4wo
dRgH7Yl1f4Dw8ajf9XbfcRDnbS1VIcxcc7k0NorDoZawjdKE0307eS6P4+Dj
n73IRdgcP/J+fhIH3aL8OzQiHpL9m7O36DhUbOJ+9pFwUTf/6uo0Do+sdMwe
EHbQC3wdziDq1Zfy2pjwlmpcn5uJg7eX7mMuwrXKVrNxCzic8NgpXEOc19XM
2idCyzi4ln7TciMsJiyvnvab6M9CMTZuwsHr1oUVrBL9/s/ORo3Ih/ItT7lD
JDo8aQhbqCPy9326v6NyHR0GpQIUtAljlCKhz5x0qF9nEyhP5H/OUAQZb6LD
+GpIYjxRL2ko0oXMQ4foh+zmC0Q9seXYF47y00Gkpawgmai3pmmqZ+pWOgQt
bvbocmKhh8pmirY76LAn3v/cKlGv2z9oF5BF6bDuimaLBlHPQ3+rr8SJ0cFP
0NPQkKj3dFCRN5GkwzyP16qJPQspNsvmfZalQ57tF3Ytoj8WuNLdw/fRoaNE
1krmHAtVGu/ce1yRDixnqxEOGxY6QeHJqTpIh0/nTJ5lEP1nO87MeqtFh2Wc
o9OT6F8JuSuul4/RAZ/lWOY0YyHalXFpWR06nBjLsHtG9P+1BWrGKwM6pCQH
zqYYsdCDdR/T4y3pcM74Faaky0KmeipOZmfokPbmQYTf/+dLVKH4Zhs6lJ3/
tbWGmD/p/Olpkefp8NQqc5sMsFCFeOTLO+50uH9nW2OIKgsFXmCzV/Wgw8Jr
sn8wMd90cm6LLnjR4f1N/vvXD7JQp/KVFI/rdFBl9WQcVyLiA7Nkh2A6rD+y
RcBNhljPfmeCbgIdZhke0jR+4vxexVute0aHH2fUhoy2EOuP82ytfU7kF4rZ
sol5ruDBFq/+ioiHLe7ZMU5iPf/xx/L5dDiurH6U/R8T+SQXRvPX02HEKmGo
coyJyleSI+sb6IA8408m/GCiFfuIEJ9mOqg93/vi0ggTRUraXe/+Qpx3SjZj
nspEKflcdnFUOjSqRUcmtjJRY80Fef4ZOpB3SGsqvyXuO3EzmXoGHf7Grq+9
k8tEpnePivkw6RCse+1EdRYTDegJCHYv0yHgfJCDVCpx33bW/n3MPgWfqpYF
7scy0Y6fO9v4dkyBpqCTU5YnE9md2NBYJzwFjBX7bU7uTJT+hlnjLToFB2QG
5gRdmUj+Ssu7LokpuO/aLeFgx0THFm8kP1aYAtFIm5qHJ5nIfUO3G5/2FGws
OfYPI+73j/vucfJ5ThHnoVy6hzyPVHq3t8ddnYJ8/rJrcp3zqMj/zWPBa1Ng
/u5+s2TbPEpvbxLecWsKJFU91Njr5tE9r/UKu0OnII9tS4hr/jwyfx9qppAw
BT+oDU0dIfNoUiMg2aBmCmzneCOU9swjh5/c59tqp6Dox9ZcffF5NBCVInnq
E+HPye9tds6jL0NV+aafp6CEZyrEjXcevQv8U3Omm3h/UuuZA8tzyL/+5s8L
E1MgsLLH4lrjHOIz9JEP3TwNF7uNfZts51BaHe8mDb5pgMLfp72s5pCyegE+
LzANAhe6Tm4xnkPGsrQsR6FpCPmd4qcBc+jBf/bix6SmQf9mgqOy+Bzi+GSw
jU1jGgyWPkWcbGagpSPibIEXpuHIunPb3zTOoiGFzp7b1YSvHVO3SZhGqVN8
AevQNOSZZB5TDp9GTjkWMtF102DxqMHuz7VphEtRb6U1TcPrmXdz7mbTaGnn
iGhz5zTIjjLtazZNIz6u6YuC49OQcHGb7eHAKaQ7sWHlLc8MzOai2RQbOip6
CVJj9jMQfpU3S541iewKVfhEV2fglS1+d6RjFAWHvwoZyp2FARkGi+dGPyqg
iT7fbMQAt+i2h9lb2tC9F+EOI8YMyC6vjijIaEUu5jPSb80YoOZ3pi1XrRWJ
VFcXnTrDAKGVfVqxDi3oQey5podODJCgePBtKGpGl1SezXPdZoDKjtvD06ca
kESYoP6GLAawQl0Tvm6uQqtq/jzUbAbInFZlMA5Xov6Zn91ZeQzgHvIp4bOr
QI/OvLfTLWbA425a1rm8crS21/xaeDUDlqk7pXfol6ChrtiX63oYUJhoYnlm
fR76ELHs3ENmQF6uhJyuZQ6KP3Je7nUfA3he6K03eP0GncxULDk+woDgrhMy
sccyUOWt9pYQnAFVoZrvbqwmowSFw7Fm0wzw3mDh3vUiEXmPpliKM4h4IlbD
L2jGI1mjK99qFxgg3el8MPzPA8RB6s14tMyAhbrvlN1HwtD30iPuDn8YMNVS
Zxrw6w6qdn+tpPyPARbyyTxJy+4ocfemxbU1Buy5/U1KTudBzf8AUTy+oQ==

     "]]}, 
   {Hue[0.37820393249936934`, 0.6, 0.6], Thickness[Large], 
    LineBox[CompressedData["
1:eJwVlnk8VO8Xxw0t2hRtltJCpdJOSvi0J4REm6yhtNkiQiJUEookS5GlUvY9
aq5l7MuMveyMQcxMvtkZfvf3z8zr/brnnuc8z/l8znM3mVrpmPPz8fEZkT//
/1cz72cQA3oqhz6Hd4ZZblSed0LcukNMDtmPF24kwltUDnwICU8TOwm2gdZt
0dBpFeY6afcHWzXBOPYff+bi5VCvakqNENPF3SG/y/ZvJCDtuz52bsMVUNfE
NM36bMGM2vW3JlsN8c9n47FlAnvwx7O+yCTYBNsVVy9XEpPHJbPxH6/FzDC3
u+3cTd8jKDohnlUYboEzvUlZlnePYa+UcvLIBkuISUq5b/Y4hXB+4y9bom8j
toYueHyeGgS7PT5e3HoPLZUBr7/+08T3ztqqQl1rXK3SWbFJRgfW2mZq14Jt
oMobPDyyVA9te56eeCl2HwLNsyLjTvoI+iBKbNG3R/TbZZ2C/EZQWx6v9DPc
AW4RjuYHjE2QyamU525wQsH6v31MK3PcMTRMfWryEMfreiiXI29gczV398Zo
Zxw+8u1k7aAl/BJEZM5vfQQ+OaGEvLp7OLU+Jmbghhs6DmtaCWVb41Fx6Ip8
3ceQCnGWSV1pC7mDsi+vBLvjtmfTnpP/7iNytZa3j5gnfn4O+Cl5wRnOVW+L
Ln73xAXfzyqOS11x0auTX0rfC5nWjHGpukdYOmrjlhvujSM9weZ7T7rDsT7I
cWjDc3SMWoRNqXqR77VlZhPPEVJEnU696o3dJ7eOepr44F7uk1J+u6dgpmXZ
rI9+gRE+SS+jrOc4H/j7luZWP4TJvy2o9PODrLrUF/ESPwyWev6aq/eHoMCd
PtYNfxBDS3+3BATAKGc2haobgOYbVksu9r3C6u0Buy8Hv8JPjXnZ+3KD4C6Y
tu25WBDe2qx3tOKFYseZV89SrYNAF5NO0WSFodbLaqClJAhe2mX1f2rDIS0g
+3W3wxtkXbu86ULWexTzYnbX1wZD4KbfyqsXo7D0X7D8Bt93CKAs3i0WE4uM
ffbBqj3vEHFAI1DwYBwMrS+M2yiG4vjqZIvSsjgksZfnFPWHwuJ8aHzG2Cfo
9D9TunUqHPT5tbGyFvEIaXt4ImP2PSZeXbmZHZqEY+uuxLTrfcB9Vir3vXQy
/lxVmC+Y8AF7t2YfK3BJhuSQ0/WEXZGwcTnz7fOeFFhZHtzC1xqJ8ZYZf9aH
VIiYJX2OPfQRuqpC2ulfMnD5SlTS3+EYeE2ubeozzAO3RXfxR/FYdC5MXvqx
LA/ehoLmF07EgiVXf1ZY/gfSrluJZwbFQs3Min+Z0E8I3VP2claIw/pDR+7c
KaWC5vHr6gLXT8gJMOjbm5SP/d+EF6wTjMf2L5bdPV40LJlzv3ZGMgnL/Wu9
tslUI8TykqHs6SSYx7FU862rIV0vayx8LwmuTr73P+ZUQ+VL0/WWH0m4FRE4
fVijBjlfgjOydiejvkVDWV+RDvaFoTz16mRsCFfdoWfJgN7ntxV2y1Jxwvuq
vVZnHbac5/QXvkjH9V3lod92/MKD6sMIGcmBddLhjQ5z3VibLtxOoVDx7ZZg
iZDPADK7T1TWKBXAcyStY+72XxTcipC4MlSEke+J1RKD/zC6VXl5wt0SXLC9
69XLHoNK9/ESH8cSKH/S89k3NQbvCFW3m09KEHr+0iaHBeNYs0qXK/WuBK7L
bzb3S47jIOV2TWhRCa5XKe+x0xqHY0uI/zOJUuxrPJgpnDSOGf+R5WZlpehI
cRDafnMCJ9WnSo7Vl2JL47VmG7sJ+C7ge7yhoxQVfA4NKY8msN51yd/fI6UY
vjejsvbNBFRubaaf31iGeUfbPjvkT8DtpHYAHpTBNE0z55HoJPgnv66Q2FKO
sdFdcneok5CoOv4Pe8pxVUe+uqV0EnJRvxrMDpejnS836njtJMzPLgxLPFeO
F+9DG0aZkyh9Z7rluEM5gmI1goYXTcFfUUzRsqQcS4oDwhV0prDOxft6lmUF
vpyrdShtmYK89vrTrXYVKN8apbuCOQVN6XQZyqMKuO43ctIcmoJbZRdb7XUF
9ps3P0yYmUL3OmXH9u8VaJa9/GZ03TS+/PjnO39pJajbRic+6U/j4Jxx5oXE
Sujs91QyZ0xD4xzhcSq7El0GSiqnmqdxPWyDlkJBJQwP+q2W6JhGgEJ7n3hj
JdzEf8h+HprGgJW+WDevEmXuRbpnFs4gvFPPxfpcFZ7errkapDiD1N0ZqqaX
qvBm4Muf6qMzKHVZtVrXpAr54cXDs6dnMCpal6BgX4VWS34tXJiB5nntjtnw
KrBynD4dvT0DSoHa8ZdDVdjoLzFy9N0MRFfEC7mNVWFiRdYV+vsZ7DZc1GLN
Vw3sVm3XjZnB1akSO91V1bjwc/ke5aQZpO8/FSuhVA2Lgr/6PkUzsIiGYPyL
alysum0gwp6By/D7+rA31fheo398eHgGgZiNfPmhGlLJO0aKxmZAtOQdtkmr
hmVgXrwaHw9iqxRvH2qpxqttzi3VIjxUespVF++oQS+XkMyU50FxeQX/K/ka
1O+TDjusyMPndyYK+kdr8IIZ7ZekwsOTRL9Irl4NIo3+FtqeIeOb+21F3Wtg
6/6oVO8yGW/q9qn7RQ1Mhc8Y3rnGwxr26tZvwTV4umGhir0xD8OUE6eOfash
9/mUZnCTjN8ZsfZWUw3M7500eP6AjM88oCHXXYPq0aNf9znz4Hm0/PHsUA0+
PrskWvmIh+T46ftnKXSkRIgcafLi4QVVpr1tDR1t0p4rTwfycHzAI3nBcTp+
/vcr3f0rD067q8tvn6EjV3Pdu7OJ5Pt2Yr10DTryor6eE0jhQXI2SSzsEhl/
5+RqjUweJkXaPPbcpUMzxO7DfoKHvZdlIoJs6WgUHfZ+VsCDRYRd1uQDOhwT
G8XpRTzUb1s8VOhBRzEz9KBcGZlfSUHvcggdE/eNSrQYPPS5e1j9iKDjhuOB
sP115HolVc83R9MxP6FvckEDD77nzX4OJdDho9l1xrOZh8LgpObzaXTM9Bzc
s+s3D1MtU/9lZtPRPuf7qqSFh5sWr7a5F9JhsaMoi97Ow4evrcdYpXQ8rKPc
RycPjX+3XVOvpuPXkcmwD108nHT+GbD6Fx3iLIWNckwenIlFXx+20yF5X8b+
Zi8PqfP1aB09dHgbbNnmx+JhQ8Dg5BcOHbVv+IIS+nn4aeL2JfUfHT98vi98
O8CD4YGVV3In6PAP3Jhi94cH3rxPgkU8Oq6EW748MchDeKNidiU/A/+66t4t
GOLhyOfqGw0LGdj9L6Yuj+TfTqZr25cysDx1+yELNtkf9bFiljADs7qZRfwc
HkTX+zhw1zAwciPO4TXJWZz1WyckGPhcdU19NZeHi0RKA98mBr4SJ46+IHns
1SmvRVsZCKqv1Bkh+c31X3IiOxlYM279WPsvD3Lyd5niexlYZ5Nc/IHk+gWU
ICl5Btr7RKS6SbZrDjohq8jAmMG2wLXDPIjEy/yTAwM35DxXHSU5xTnvo/JJ
BpJ2lsTok6x9Tlvn9FkGqGZvjt0i+a8kk6KlycAK07D+2yT7/32QcukCAy8E
X741Jnl3wRIT48sMvIvcoqFKcnXghxWWBgzwN22bJ0XyXfMDhI0pud8fBj/+
kvUsVSixeniDgY6qEPsUkr8K6m94coeB26yUHeYkq/3mVL+wYeCK1vPfS0ge
+OrxKMiBAbOcBU9iyf0/c12zO8KZgfUQ27yf5G1a8W2xjxmoMs/KSCHPs3ij
ystELwaOS3YckSbZ4j+GUpYPA3G0N6nPyX7MLzIfovozoL61f2032a+YN5Nh
pUEM2Dj+sdxF8skbL9UZ7xgQvZ4dd5vst8fijPjuaAa4FWq/CFIfm1pVrw5+
ZuBO6e2KJlI/RELropEEBqz0haO6+niY1Z5nOT+bgYCto9wqUn8Rm0NEhX4w
kGxVYZ5M6lNpZGfpmgIGom0vZnv1kHp9e2GbTCUD58qXsgVJPU+0fexV62Tg
W8hOSxvSD2+TDr650MuA0ZWYkD+/eDjoXn7y2h8G9HeHhFwk/WQv/V/03REG
Xu81nbeQ9Nt/t46ZBgjWQlFS0aKsiodXSg3C75bVoitHeKS5gvS7kGV+lEgt
DNYvUm0h/XwvJWBj2rpa3K3IUs6m8TA43tHesLcWn/dcO1aXx0PmDf8CtYO1
EBb+lXrqOw+Pm1TiqEdq4a017P85i5xnWRF3v5yuxb1vZ54hlZxPDvozLtdq
MaQofiE1jodlrMWdwya1aGJm8X+O5qFZ73uhxY1abMyLXe0XSa4vL/ZC27YW
ltuFTXeG8hA60igm/awWjZk0t3hfHszNvHkhL2vxYlZRTuo5WX+9fNeywFqY
ntE47UvOx+K0oM/jEbVI/7TLfQ85P4dtdRQq0mrxxS7s1Pd7PKgOV+nadtSi
91BkaacGDytNXA/1MWvhbFK5y0iVh3a67Lprf2oRYTY8VHWCh/vJL3pOjdYi
N+iTkwt5P3ywOmsrtqQOvmKBteu3kX5kF/kTB+ugyF9ca8ObQb7BfTt5pTo8
fbx2ZnR8Br5VUpfij9XBbs9a71v/zWBzgodkkEYdDucoUKX7ZnDuztGEG6Z1
YL9ubWurmUHMn9xyIb86LLhqk3CdvC/vvdWSuhVYB/NxStrSkBkonOxxpoXU
4fLOl9/jXs2gPGLxLpfoOlxQG6pP8JzBsPYV/4HsOvwz0DYUvDUDlayxC0U9
ZD2pg5J/983g15N9bU6H65Gm1bbpbeY0Pu6lyTeo1GNvWnny9cRp3Gm77Lf3
ZD0aat0fbYqbBp/CY/Rp1kM37Vqn9ZtpbP1THaVrVo8jnxQXGNlN4772HYs9
/vVg6xQoOctOQ2j9Jy6TWQ+bI+/VSwOncPTik7aegXr8OfQt2spnCrb+RhXd
nHp8dwp0Xuw+hUaK6KfOiXo82l6dKX5vCu9ZzwxblzZgi0hwUuyZKexOvlVV
J9eAtrA285sTkzh3YvfXAs8GtNKOSWdrT8LXMsMiUroRgSLzirXbxvHnaMKz
JTsaYXnPybivZhxnRGPjHfY0olq/9rp1wTj4S4I4GoqN2LO2KO/Kp3E4Sd93
mNRqxHX1j4Y+1uOw6NjvrePcCImlTzJv8o3jmF5y7LzaRuS6mYmdXjOGsaNf
mTfdmuBqd249/7oR3N9flL6msBl+70ZFR9q46EgsesHY9BvE9J9TkxG9mL5V
H7b9YQu6jlZvNzrXgJ3HFs5a5LaCdeCc0fuf16lnPm66kbK2HY/S/XKulTdR
H6w32tx1rQOx28r6CopYVNkhvSVN5zoRYKGr2ef8l7py+V/BoqlOqMseocqE
jFLnUcOSpMK78GtqKeNj3ST1t9JJ1Y73XZg/l6L8mTlJTfo+1Bka1QXJLI2P
UaOT1KuZyiIin7qwdlhsn/XaKWpKQsd9vtQuCCHLwfnqFNUwXEqxraQLRb94
Z060T1Gznb4WBf/rgt5cYolh6zTVb/yCgc5YF5bfXS5TNzhNvW4/M7pssgse
YxUuitPTVCGbc9u8ZrtgMjHB7BKboZrf5D63W9QNtVHVttcXZ6gilw5oaW/o
xt3GwWGHqhnqbbncX4vUu6H0Q9JjdQqPStAbvmWf68am8gfioT951NV3/7rd
0O6GX8e+U8KVPCo1dstWml43fgYvONzE4lFXrg2wfWzcjUiFOA4hMUvNnTRb
Mm7fjbg1ktoLnsxSl1CFlJmR3bAKDvXfe2aOaqy/fUVgdDdy0zNXzOjMUTPG
T/Qci+vG4+HauzmGc1SjvU7PP3zthtltm0+LHOaoaVE9DVczu1FtPvR8LnqO
qu+ZfY9R0Y0NXvIZSrw56tezppE/x7oh8Aotd3X4CDnZTSqY7MbOBr2uNZf4
iJ9CnS3U6W7MvjR+nKPPR9DrDNbm8/Wg69F3j0EzPmLE4Ipf4eIePBwVus3/
gI9QsdV2KZXsgcHbYsupd2R8qMrl2tM9kO+5Xi3YwkdcceWN6pztAePiUkp3
Ox/RbZQXWKfeg34Hg/fp3WQ+acWaeu0eBGuuWXX6Dx8hliR/uulqD7aG3SkX
m+QjTAtl5Vrv9SDvvoXpk9UUYnRQfAUruAdlUsTIpbMUIiiiryT8HZkvSvpc
hAaFOKCV7nYhvAceoSfzf2tRCJvUc1wiinwelzOhfJFCsB0fV4d968Fz0Yjn
saYUonc+y/d8AcmLn9+ZcKIQnlmpJwVpPdD4J7GS6kIhpCzdZn6W9CBJNrHE
1Y1CmFSK3pGt6kFPyrLLPZ4Uou21uvrC5h5QDjgJaPlTiIYNKYt+sHswnHhI
4PNHCnGf4Zpv97eH/L6QOsaNoRArn6g57fjXA44LPHd/ohDnWT39byd6sD/o
ieibrxSi6tuaUlsBJhKOx1nxp1OIO4Y9btsXMOHYYVi1KJNCLFmRrNApyETM
u7xdi7IpxFm7s5/OCZFc7jbKzaUQtMMu3jJiTJRRWGEGhRTCbFBVpUOCib6v
5ou20CiEQMTqsTeSZH4zJUdmMYU4xpdkLiDNhKJF3BX1cgrxo7jrZPtuJr7T
vXfk0CnENcfEmaB9TKhrZEYdraUQ09ud09XlmAjvlRUn6ijE4ZerpHMOM7Fo
7NTyT40UIuPCGYGgk0zoZIUtW91KIZ6pZ3utO8PENtk5wTNtZL4T2wVjzzLR
tqd3mW07heA/sGRZpiYTDb9m5dI6KYSmSM2a5stMoD92totJITYtPvrOWJ+J
80FHNVp7KcQIJUViwICJs33qsTUsChE2/HrjlCkTFNV/D6P6KUQ//eKOdXeZ
uJa8/OfwIIXILS35FmPFxNfuD9bpQxTCjzi0Z5ctExXGj/dZsymEXLK4nMoD
JgJzBzpqORRi4ecXmcVOTEgWXKp24lKI3x9mDmm5MMG90Vkp9pdCuPu3Kxu7
MxF5SX7m5DCF0H2qRfQ/YWIdd95uOskybsRxG28m4ohRG93/yPNy2EebfMbE
i9bXpQySq+99POPxgolHa27vV/1HIaIsVpYv9mNCf75wchbJ9w09NQIDmGDy
SxzfOEIhxDUtzse8YWLOoCS2lWT2qaY62RAmNI58cdw7SiEIZdWLGaFMyG69
Z+pKcqB8TrNyBBN/j1OvF5JssWuHfvEHJnY8OOvKP0b2a0tYm+ZHJswtfyQp
krx0/VLjphgmqh/XTN8muXOVa7fRJyYGmQqmwSSnLeWY9X9hQsU1tSeHZO95
Rn3W35iY+Dv7uJHkqzM1lpOJpB7iOhWHSN41cnTIPYWJz4WSItMk8w2l3Fuc
zoSSrOdCgXEKUdezefh1JhPpHh0S80iOawm0k8gh+7tNQHuWjHeqmzcWnctE
vQItZphkjQp7R9mfTIj7z4m3k7yhkDWVTjCRWeSfVkjyf98vuSoXMhEge8nq
I8nFqaVzNBoTtxsVtR6S/C7+sIdmKRP+L6V01Um+8zF+XlM5E2eaeY9Wk3w0
VOKpURUTHuM59GbyvFa+9l3UX8OEHktd7Q3JrOe8F9a1TLy88nVIjeQc93tC
k/WkX6h12ZNkP3ydOgLcm5jYp1KUFEWykY32ysW/mQjLcKs9TvJ+y/w3r1uZ
sMtburmN7O98k/2iEh1MCEbZf7AmuflydGh0FxNPbpaq80h9uKl6fUhnMdFz
wkGen+Raw5vX9g4wEb/s0H5XUl/S9upi3waZSG6SjBsm9VgeJRwY/ZeJrfSI
8WJSr+uyR7Q2/GNiLNeoaRvJVtVNS8NGmSjcFLrfg9T3qukI79dT5P7Nho2k
SP1bCLufEOIx0ap+Qe4W6Y+cbWZ8PnNMiBZkPYsn/WOou+Oh+7xeKPQpZK79
QyE+J2TetRbqxYRe5KkY0n9The92DK3oxRcJnyNZpD/P/Xbpu7GyF8Mdu3IL
Sf/+t+CEiZFoL1J8qdfKuimEknGNrubmXlz69/LcU9L/fg6pwuXSvRiyTs2w
IOdDl++b6lPbemEvuikV5PzwztFXVZbtxQIh87vtvygEXaT/iOzBXggNz6yq
qacQm7dXTHw61IscL/ppZ3Ie2SMxQ+pIL5L79UQ3kfNK/M79PRJHe9EUYhl5
pYachzQ+qcVqvXi2RdTqURmFmHggtrjfoBfcgjCzkDwKoeY3U3zduBd8z6fV
Or5TiIiYjicdpr14rGjwZ0MOhTjOiOU13ujFw0ucq34ZZP937B+m2fRire+D
AxsTyXnVptYc7dULv2W3NMXeUwifgiNJes96MWKqpVYdRp7HJ1nvhS96IXjj
r4DLOwpRYCskdzugF1Gb13/PDyLrE6z13x/Wi/kSmx+M+ZDnfeDKmfxkcj+G
T28tdKAQK8TUJO3SerElhnPvrB3pl1nFUenMXrS3jV7xtCbrL1sX/Sy3FzGX
sznMWxRi0LBzVqu4F1736EbnjSjEbp8bme0tvRgcy63/dZpCBFtdfhnQ3osN
YkKvvp6gEHO6Z82Od5Hrq6qbOBwlz3vDTpE4Vi8WydicHT9M3qcZnLt3h3uR
d6/QN06WQqR32m+ZWcDC1u7jAT9WUAhFBc8gsf0sxAjv4PTU8RGNpmIHI+RY
KHQtvCVF5yPs/BKbNiqw0LPviaJBJR/xrbdZXEaJBf/cunZqER8hGbQr6uBp
FgYWTPocy+Aj+IcbEy5cZWEFcZ536A0fURG/veSlBwvK6XtUl53nI242/rwp
7MVC+7gVM1yDj5jHr7v4zVMWLNW85bapkt8nVx6dC/dlQcd6BW8n+IjkhXV1
8cEs/PdF0WdMlo8IMnPpLIln4VCgst3rhXzENcmaSUodCz94Q0umbOaoQmK2
GbkNLAik+O6bsJij5q9cbWPfzEKtcnb3H/056tZF+v39bSzccJlbknhqjsr9
x2qs6WdBMXTO6ZvoHNWjbDYtYpaFVFfBX6/yZqmf7+++d3hHH4odA/OPzvCo
o+Uvu63d+5D9c+zxRYMZasg8w/s60v1gx15fKDc1QV3ZnvGRv6gfKSsLjWPO
j1Lnu7rrHrYZwKHhPCZXepjK98UwqVvkD07f+zaWZ/OHahua9FKk4A82B9Q2
HD7ZTf2x85mg8L1BKJU67LZpb6AerF9b9dp6EBv457upfG2gJrt8erXKbhB6
54/3yjg2UD9WFYuLOg6i5HKCk55IA/WZ1fxdGzwGcdms8J6waj1VJ83j/K7g
QXhMn8rYl1NL7VN0DT37k4z3ONXwn3kN1bhnqVEFMYip1vgJnYU11N8vwqXO
FQ4iU2Kp7ejnampla+437dJBRDYZmPKzq6ipj6Z+XqodRJz1qPAJx0qqS8GD
HnPWIAzW8g2/CSujCqvZynoIDcHe+anHzrVF1Mj85UsUhYewgd4SsPJ7IXXv
4YSB4ZVDsN6f561jUEjVlOmPMxEbgnflfa3GmAKqzwLDTcekhzDvQuUFnkI+
dV7h2TV8ikNI1lNKXNGeRx07sonvkfkQFulI03knUqleaT/b5W8O4f1K0/Kw
yBTqqp3XfrBvDUG77NriSV4y9YD424cG1kPY8bCdzrJPolpPLB1Vch6C25Uf
aYPO8dSB9ImB6YAhOKxdb3+CEUlt3UWvc8obwhZfOYmMx+b4MCjsKkAdgj87
LruvwRamXy5sfZk/BKXo6kcHOh9iQLrJMbJ4CL+6Q1ofWHphTKJ9fQl9CAeC
X0as8X0N4cVDN1b1DoH132BgfnMkTrMWziQuY2NEuaS4nZWIRTFnYw+tYOOC
+zK2XG4SKk1eaBaIsJFQNvR9wZlkXGgTimpYy0ZL69j7qLUpMK5ffXpmExu2
7GumzUmpcC6QCjh7kA1ORXPdwrwMJL+HNNOQjRsWl48+OJaLfw//LTQxYUN4
nZhj5KtcHLz0abDtOhunbxVWLOrKRe7y5WlNN9kwyt3+O8w1D8WP249V2LKh
77JRLivhB1pMXY1SvdkIXq5s/m6ACknsO7HnORtecgpl9WsIGEv0bv32gg15
mvAe50MEWHUanNgANp4rTT3++JDA35PrXN+FshHanWx/eYrAgm25oY8T2di1
ly6j0ZePswJWj3jJbKhx8n13CRTAt2Oz6cM0NiyL/8v7KlkA4RCf7fez2ahx
CNv+SLcA6xZdzb5RwMa88D2NjrkFMOpdFs4sYuOoYrCSSn0BPubnu5mUsMEo
eFb0YagA2x5uP6NfycZ5c9X+g+sKsX9wokGzkQ3tKwqX0hwKYV/yLaeimQ3Z
wAuRl3wLkR1tHKHawsYSeqrV26hCKF8rNTveycZ2xr7LA+WFUK16+0/+Dxth
ChGG9aJF8Pmi3pQ6xIbkg+csnZ1FqPKa+76Hy8aAHZ/lHeUi6Kjc8Ng+wsaY
zZV11sZFMEg6KLx+lo2bNos8sqOKkHl1QojOx4HB5fVNeqlFEFrwfekTAQ4U
184lOhcUgbimLDggyAHl1o73xl1FEBWcmx++hAM1wZpXG/4WwSaNENAS4uAY
bZywmC3C5sUn59JXcuCtoef3UoyGhxnzeTfWcCBcHrfk7lYa6oxLpsTFOFjP
y/Co30+DZ5ba2GNJDnzSVm08okZDq+nSkQObOLgYckpGUY8GeaHqYZYUB3bS
42zCiAaW2Xm2xnYOHpRZO5nb0aCyYuXg3E4OItq7pJ640PA2t74/dTcHBcQW
PRkvGv5aBLPM93EgNs+aT+8lDaoil5michz0RJ4TW/SGhsgfYt0VBzmwWBQb
eDKchsmbLR2PDnMQjMzrC6NpiKcatjBVOFh7p7hpcxIN/Lc3/np7jAO30fkW
zuk06K/pblQ7yUFN+GNJ/Rwa0vOj63mnOdC1URjP/0HD0rvmtclnOdgb4MVK
yqfBTHQb/boGB5OHNvZI02j4UdhftUaLg4Nda5gbS2lYYxVfUXaegzVn6zqj
y2m4J36nzEWXA18/MUZCJQ3FtF0ley6R+Y4JJqhU07DBhlvUfYWD+/aidy/W
0OC4LqXgzTUOyprWLx4lmV5iS6gacfCKetpVnE6DjJ3cz2kTDtqPrSOKyefu
kmO5iWYc/DvhUTlJ5vtdlpVjcoODgJdZ4WlVNOy3d8padYsD6Ra9nSMVNPhs
PJJRcoeDz9zzdrllNPRUzKQ+tOLANG7aalEJDUoPfibvsuXg/LzZde2FNARt
fpzYeZ+DqrNGDocJGthVx74FPuCAT8DdUSKPhtNOAvGnH3LwRmzjRs8sGt5L
0z5NunDQyNS9ZZdKw3iNd+w3N7I/ioG6A99o0HZWjTbyIPX1KL61N46Gz1sX
R4l4cbA7ncZnGUkDpbbiPe0pB92+t6m272i46voy3NGHg9uOj9YKvKYhVUYr
dOdLDjTTlQTFfWhYUr8ipN2fgwX6R19ludOQtyMo8OQbDpbMqt/xtqJhVaPe
q/G3HHQt0KVmmdNw132tf3woB9bJUd+s9GmQbA71WRHJgRwErtmdouGxV5R7
azwHJjKF6uUraVj2YKGaSQIHGvbxqusX0hB6864IK4mDXYVfohdMFSFN/XA0
N53UK2uBR3h7EVjCjCJ+Kgep0Qd1Gj4WwVZA4eXTfA4YqpykltdFmBsJ11ta
xEGflpC/owfpx+abrNVlHKin3ey3J/2t/p5PcHsdB3rl4llua4vQ7G9BT2jg
wOp9XpXivCKYuVeG7G/mYHN/aqjv30I8MgvZrtTGAUdH+8FISSGSd+xV1+rj
4L13okyjTSHWZBn52U9z0FJ/hDqQUoDoz7SLUzwO6sWTPFeFFGBv6M4Nbnxc
iMlc+dDoWoCzruNJz+Zz0cKXN3+TagFcTvgzwpZzEf1gDx+9KR/dNT9XFUhx
4fQtsvNjDwGrfOm201u5uJgocelsIYGZVJ/YChkuJC8nCD77SGB18MWDDbu4
MJs8bVxtRED1GudSvwIXiQsC/ChfqUjoXx8mpMHFW7EyCTH9H3gW4WXcrsnF
mlJj11vzfsBMh70l8TwX8Z8m5PK+5WFdXl7yuUtcfG9UCzSZyoWPv36xrykX
wW3V2d+ff8fNg++GFztxca1FsLPdMQsnB/myfjtzcaIpTf6CSBY2RN50iX/E
RXLS1qmK+Ew0Ljq0UM2TCx95xr+K5gycamtc99yfC8OMZ7+WbU/HZs9Vqgvj
uPCM1/zRFpqM2UMuy5o+c/Frsf6yH9LJ+MXuqY37yoX528on6vpJCLiUZnA6
hazvosjN7U4JmNuhY+eVx4UtY3ULM+ILWhn+7wXquKjOU1w//SIK2d7j1+sa
uJjUT9U+4h6JoCNG26Obudg4pdt1YOA91GN3px9v52L8dpX87a+h+O5YVeY+
QO7vp9q/Sf7XCN4l739+iIu4DwLvJD74w6Y7XHcTlwsjxw0NjXtfQEbjTgcx
woWv7vtLq1U8MI9SHxMwzoXUpFKFip0zOjOO3DKe4qJKeJHYG31b5N2K3rOX
x4VrlMk8y1umeLthyejcHBd3rE8r/5l1V/kfCvSJbw==
     "]]}, 
   {Hue[0.6142719099991583, 0.6, 0.6], Thickness[Large], 
    LineBox[CompressedData["
1:eJwVl3k8VO8XxzEtSlESUUJKkq0S3xSfypJoISQprYiSJYoWWii7SJJkpySR
CCVzs+9r9m2QZcYylrGMit/9/cPr/ZrnPuecz/k857lX4rLdKQsuDg6O2+Sf
///XtRiqI+jG6hNaDwPWimxTW6IhYt8trISSTydLdxE96nuiwiK+CGtCYcUW
gXt5nPi9aeujO1InsOrn/Y8vWOugV9Wc/lbYCDeOycy6fJfAVj/RhEUxU9zr
EQqPTN+Jv7pXXl2SModkMNXtJXMPGB6/Ci+FXsLK086C2uL7YXJ19kew8FWo
d9x56Jt5CIUaIlkFEZZA4C7vLittKEqqpbHErLE6zZkyk6iHCK6LSdviriPP
i4g2ztEHd+/j2NNSN3H60T+P6/2G+Earryowsse+xA2V701MYK9/VfdcqAOk
9WROsPnM0KnwTMNf2AlL1x9ccijpEkKiNhDbzJwhbi2mpzN2Bbp8Hw7kRdyG
y6l1ngqnLfF1rHIvU8wVm0yXfnQMuI4b5ubpzy7dxVceZmq2+U1sqWbKi8fd
Q+EGR5vW9/YISOGXNpByg6XhtY0VNbegJRofT7dyR0zAzQNKD5zhVhy+5qfR
Q9jQxW7R9t+BkrKsv2noI2x68Og6Y+QuoteffOoj7IGL/gds/QMe4V7Vq8LT
3zyQ37Yxh530GKc9aVySZp6QWKa3MFb1BKumHdy/RzzFOzVNPyuFp3D5FeIy
IuaN5yvp3kpnfGHo1/k1m/CG9AqVT35ZfpDXlJr2uOSDI4ckJQ3FAvD7S5aD
aJwverSJp7wez2Hwos3mhFQAnCdkCPGoF5DVk0wSKQnA6mvv3Y3sQ8BNuTE4
YBUIU+MRns9aL3EhZ+Ez1eg59Lw7+Tr/hmL9jufyZ0KD0N+c+6wiIhyPuL9s
9xYOwYrNcmEJ+jGQORLklU7uI9uVeCC7Pwb1nnb09pIQJCZ0dux4EIutFNlk
+dsvIf0kNWnrlzgU/4uX/1UfisvTwXyn9idi1VToXjG/1zgQ+WG3/vdkZO5y
DtXpew2dMG/uZOePMLc3nHVQDQfD1mbVid0pSB3lyykcCoeVWm6gZfonnBry
OmCjFYGOFZd12gfSENZ5VyNzIRLHU1L1clQycGiTaXyXcRScW1ip5YkZYJxV
WcqdEgWBinP3pwQzsXnE9UqKXDT0l4XNPJvPhJ218jaOjmg8yZOtjK/OAv/V
1PcJ/8XCcnbXyEjSd5wxjUkdn4jH4PFjb/9blw9mu9HKWJEEuJ4RrAs6nY+n
5twWhhoJyPkyr3omPB9frtiJfA1JACU8XXinZAF4b6p53lNJhP7s89LOfYUo
etx6dtmDd1j9pqhs471i7P64dtkm7g9oMmQ0r9KtAM/io3NHNqdim/FWviCO
BoRZm5jLaqciTyz84Em5Bmz9JXtx7c1U1K/VO6F0tgHqSc1X2n+kwuXOagHH
jAbkJIVmZsmnIWOXcWumzS+MGo7k6lWnYbXlEXOZrkYYv39VcWt1Ona8MToq
0d2CbQZjQwW+GRBtctfN3N+NO9X7EMbKgdL+U5p/+QYhlLG2i5OTis3VcmdW
7xnD116NypoD+UgsHZX9+mMSTTqPJvyP5OPhG8svwhWTmPmUJ3jsVD5e7t5+
w6llEir39l0qtcpHl4GBFN/UJLLXKcwQwflQ8Pr3pnn7FHI0RcTTh/Jx61la
dFzgFHITx2+FhBRAnzP5Y6wJCx2r5F+fiirAQVqtfskVFv46Xs9b86EAx/47
WEqzY0ENA9z+1AK8fNo8OPGMhbzm9reejAIEjFqO+GSxQOUuKbl9sBCyImWa
XgLTyLd5u9F0pBD3BS7fUiqaxhf69evmM4UQDqsvE6mZRvw11e9XOIrwuqvC
bbplGp6WzWdvChRByn2q6tnINHQurw1/cqAIQgY+Tzevm0GNqeeGT35FGP22
4xO3+Qw6dGzXU+SKYV9ncd2JMYOqkv0W3CrFiOhLWaicnEGe9srM1YeKcbH4
msGGPzOI0nxntMG4GBLMPUue8czi8sGeF3IPisG3wl+3YOcshlSM+c9UFeNz
eWH7R+tZTEup8aXYlkD/7RvvybZZqPceLvFxKYFNnjW3SO8snr7Vcb/2pASx
c97qivRZCAoYMSVfl0BJf+619OwslDmv14QXlmDx5EDQEf45uLSHBXptLEWW
R0Pvbq05EKGRRyylSqG4pywQenPgPhW/qLGrFD+XT4urGszhdWmq3aJ2KXzX
vP0xd24O3zKLT952LEVF4/iT8Ftz+BvI4rtaVormkwejLN7OQVNvvuTQr1JE
8G+o1Y2dg98yjodi3aVQyD+wVfTdHEQf8Iy3sUqxmPiT2zttDuo2W2oNxMvw
5cWeZ/fz5/B0m7SXws4yzAVdCy0rnkM1Te7gauUy8EU/saNUzMHcZF9aqV4Z
Euq2Gqk1zMFdU/857pTBesf53NLeOZQuGuuIPi7DcJon8+7AHNZ8N+P441cG
Syy8EmbMIXqXlf3X2DKwOpeJb5og69/spi9XXQb5oaZJpQWy/rYn3DytZZCN
NPl6lJMN/ZfexFBfGQ7a3InSXcJGD89LxXh2GRoX9uSuWckGFzt5zcZt5Xh7
7ZxQigAbG6sOT0GhHH9jNjXOCbKhFNPaeHVfOcyWtC3ZI8yGxdHlbz4dL8cj
8/VvHUXZcBONfNBgUo449c2zd8TYCJ1Qujh3qRyeW1UTr0uwUfr68rbDt8tR
KmQzKLqNjG/LXm71sBx12h7XO6TYmD/0nOHrUw6f8h/S3tJsyDJ+pDZGlqPS
VGtPyk42tPKMguffl2PD3fW3tsqxYR487CT2pRzlkivqvOXZCFQVVrUuKYfc
vZVVm3ex8Z43bVNAXTlG+gTNT+xmg+jVXkxvL8fMac2p63vYaP3a2dPcXw4F
g31P7yixMenjVPiXWQ6xxeQ1DnvZ4LnA805ivhy7Dl/2O6PMxtY9sd7aSyrQ
kblrSk6FDbXl+25c561Absai+gTJp9trTjzfUIHI8z+tY/9jwy7VclfmlgoE
hN2zPbSPDa8n/9a1yVYge6OGZjXJ0SYhMwvKFThvf2BIR5WNnJ07WyUPVeDd
7BvDzyTXL/78rqNXASdmwsPl+9lgNJyJtDWuwIY7n6z1SKa8Zz4MvlCBFSHC
Kx6QvOn+0ytZ1hV4qmxmGkHyXn1R7Y5bFVBUnzVKIvnE1gxpTrcKeLk9mokj
2XJOl0fKqwKuBW//8yPZvbJnVDe4At8bF9dfJjks2qXWLqICKu7XvbeRnObE
9yUksQJmgRn+TWR+ZTqJL3PSyPy780SdSe7dpObS9a0CzqfdlDhJ/jPecJZS
VIFlYzOVd8l6BYps1KRrKkBkHujtIfWRe80pfry1AgnNt+yUSda2DeNy7KuA
UEStnQup54VDCv2hoxUotH3V9Y7U32V9ccn32Qo879yUVUj2J+nHlN/SVZVQ
/LzsRTHZz/wgHzsZwUqobLDwTyb73WYhceqkOMk7b3c/IP2wivek0GulSghn
nHZlkH7Z1tvP/qFeCW1WTLYH6Sf1r/c7enUqwXH/lfoqWTbszT/EyJ6vxAhj
J28f6Ufv3Yc8DKwqEdFQenL3djZil7VY3naohP3qVaU3ST83fFoqS3hWgntB
6Xgq6XflxYtfDT9VgnnG4fAJ8rwcO0481squxO8nm+7yC7Fx5Y3YSZX8Sohu
X+JVSJ635ypdgyJNldizdrnIJB8bdDsz4d5/lVC4HZPvuJQNjrxv/Q3Lq7BF
188kgYsNwVUi6UVrqyAjtT+5aHEOh9+36CZtq8LQyFGPGvYcImjG9+2PV+Fa
lubX5aNzSJfP1LlsUgU5hQX2Ozo5X+4LrDe6VIXC2/VDe8n5Mr2hIUXFuQon
DghTRbrncMJAv3shogqHfV8LcNfOgTNf97D/SBUeND8us/s0hw1rPvC6z1RB
PO5nKz1pDvLmK9rtOapBEdVyNUiYw9n5kltGAtUo2a9+if5mDhm7tRI2HqiG
Hd3w636vOVjGgfuDbzXEfK6nW5+fw/2JyF9vXlbj+sBarVyTObzAQrR/VDWG
ubzO/iXnN9Geu8/hSzU+Z/as0tKeg7CA6vX/2qsxqhjpOyc3h0oPpepimRqc
PGO3bIg9C1W+Cq6gvTVQfBC8/8HkLN6/vqRidrAGUYfk0v4yZvHkU0A007gG
SX4uZt/ayfUtQ44bHtUgZKS1RTGXXL/zrZBNcw0iUy3f87iS99HXPceUemtg
r5tqp2U/C4+D5Q8XRmrQbtoYYG01i7QPf5yOctbi3o/cTnvjWfhSpbs6BWvh
X5vOeV9xFofpj9OWHa7FnUspRU97Z+AqX11+/UgtONLXa1m1zCDtlnB/7bFa
DAdQdHdXz2DzQqrwG5NaPDmdN/ksZwZs/s7HCra1eDxlL707kFx/QMX4TFgt
bIarUz2VZyD2fJidNFYL5QTjuqPO08i75J6UPlWLzpLzc+bW0zDfs870+1wt
/nNSmLx0fhoRTarZlVx1oPw44LhLaxobRH1uMwXr8NUof8/sumnwf5CeUkId
jj+/ZLvlEwtLCy1GqIF14MwwsxZvnEL8S/ab0pA6GGbbFswXT0HTyl+v7nUd
mhQpnwqyp/B4ZeaH3rg6vHe4pi8RMYUF/SXWS7PrcPWrQuCqy1OY64zt16XV
oeUIz4kkxiSGZ7u7GhXrYfZgVfQu5gS+WgXm6yrXIyzwwLqAzgk8bFZPpO6v
R/TyBYXmigkIZr21TdKux2ODimCF9xM4fNvs7/1z9VC8NbZv/sIEwllNwlu9
yN85OxQ5K8ehM1Fl5NhdD/8bex8ayzIRz/hezhvQgPPqF7Vjh+m4+eqkpM2L
BvRsUD2oTdChotl3ryisAS4BL20qQugof7tS7n5cA/zUzpS4qtExoW8aSM9u
wH2HC/3W/kNQz5oxLOxrQCuf3X/p2wfR+mRXp+u+X6j67mbxQ/s3eEXfMX//
/oXq06Ypjdc74GedaRm9tQlduoXZ73jTMHMw+fc192Yc/cKfvEG4heq0uzBD
sKAFxE7DHf1n6dTuT4W+dRJteLQpVEMhcYr6x+bXmx1321H498UtlWE21Uj5
fJC8ezvedjZki3HOUz9yDjzd49GOl3MFYX8E56nmYbMOagHtGO5dNeJ/eJ76
s0hExyC2HffzfPZtfD1PfSZ+ieVS3g7mFRmTCo0/1HVNo8dKRTrQI0NLEQn8
S70Rc/twlVgHPDL8A03j/lILbyyq1G/tgEXzLxWfrL/U25S1kh3yHdjVuj2e
6P5LbVVUYjMPd2DP2RaLhwr/qFG+dxOEyLqlzlx0Q+U/6s5Dyxcsv3dgmZt9
yvn5BWrah+TtLUQHDpqul9HmXaQqCegbHC3ugP5W/emNWxapBwbD4nfWdeDb
7SHlwKOL1BP+MnrjAx0wT81ICH61SHVsPR7myt8JnTfcqYtCHESOQ8huP+tO
dAnWOnMlcRDq7f+d+3ezE+kXa/1zP3MQ+ZqdnjedOvFM+q/pjW8cRLmQVKu+
eycaHWMmP1ZwEB0/stzXv+xEUlbopc+jHATHyvbKSKITdZZ3lH8qcBJHYiWs
Pgt1odyKJ6g6mZNQCuOVKt7UhbNnD3FnpnMS4gF/frdJdEGCfC8MyuEk5l0b
Ly+R7YLB0Ct1+RJOIsXA29zkYBd8N+StmezhJNZxThj/u9YFlbPTTjHruYju
C4Tm0W9dsPCJ5H/ozEVUGqdQzKldiH8xU7X9HheRrRee71jYhYqGy8eLHnIR
wSpOByOqu0C/1ejV48tFaPFJH2D2kvHWmT6NjOEikvMC97zk6YbrP5fT2yu4
iDuiF7b0nOvGoTAi3UyQQlzWVzzdcakbP7n/7owSoRAnnnD6NFt2I813rUXL
ZgohRY+bqLLvxkbBdB7Z7RTiV8YQ8c2jG5l1GV0mKhRC8ZjjhZCP3XhWLmFI
NaYQm9w1XgR+7obwTTW/86YUYnm6QInP1244DEucmzhHIbqEsuQfEd0Q6OPK
Yl+hEH59f/7Z/upGxOiEsL8DhRi66xlx5G83vhEGX3J9KATj9M5rMRw0MIUC
T33xpxAXx9dxiVFouPadbhn1nELoSfbvFeGmIds5aNAklEJs8X72dg0/DWtE
vQV1YyhE2BZ7lQABGuglZ2VH4igEb+6ZOh4hGpRxOOlJIoVgj+1YunwTDVJd
nXXByRSi1qjyxr9tNLyP+5G+K5NCaI9lLLsnTXKOW6NbFoXIffY2ek6GBtcw
RcOfORTi3bebv6YUyP1H+0yl8iiEm/jaA4x9NPTpRNrZF1MIVg670eoADRdm
ijSvlVIIG8Neu351GrJ2RtoYl1MI46df4mgaNPhb1MbxVVMI2RFjnubjNDg1
zaZQGilEttExqo8+DUsm24ZKmyiExo/Dt9QNabi5P/rpoxYKcTZAoT3hDA1y
TvsHKtspxMDMtkBTMxoUaLSEC50UwvHCJo3V5jQYOC8v7O+iEN6KK5KdrtDA
4zwsVNlDIQReL5pLW9Kwa0mVnmwfhYjmnOHvuEbDzuUube6/yfj1vXc1btJQ
rDy4hHOQQmjub5WftafhwdGtUfJDFKImrqb3wy0a3N5uCjpJJ+M75eryu9Jg
5ZLgen2YjN+ZvlB0j4aG7xQ7yxEK8U8rKd3VjYbn37JzjEYphNenKEu5RzSs
b9MyVR4j8xEKFel5QoOwqMpxHiaZj7tfdchTGliSki8bSJYdevxYx5uGWdlk
xcBxUo8cO0ZqAA06p49eok1QiGoJy8grQTQE8jRPO0+S+vicOyUUQoPhKLNj
geSByVPLKkJpuG1oI3hviszP7Og3t9c08O3dHssgeaEAN3dHkPpL8z88ziL1
klXeMhBJg/EVibQ4kgVeyja9jqEh6Zvx7lGSo/9t8TkeT0N67LdlMtNkfpbC
6pzvaOiVt5Q7S3JWNd9kRhIN+WtdEx+QfFhlWeK1jzTURGx2DCG5Ouqv6aZU
GtQJv+eRJJ/lnlpd+5mGppVLuCJIHrCn/3ySQUNIUk2xH8mOrd3OKlmk/lZa
3fYk/zvUtGM4h4btKUV6OiR7fajsjMwl63v4fdU6kgXWFQSdotIg7hGxpe7/
+d7L0VqWT4P27JD/Y5Jlf6eycwppGMvYeXI7ydnHEsnPdRo898xbUkk9NDIj
LkmU06DrwWrUJbla9MX6xkoack8VRZT9X9+n3mVeNeTzfRK5+/+v75j7gwP1
NPxeX64UQ/bDweT2rvFfNLSaey2fJ/tF3ij9cc008rtj115tkr2kr7w2aaNh
HyUxz5Psr0CQ6XGeTjKeWk18NumHKPZJTmo3DVNHY4e6SL98LT9gLdVPQ+Gf
4/5LST8d3rNHtG2QhoPiWuNLSb9VvdlR58+ggUNAMpNN+vH3DcF900zSP9Ga
1t9J/67jG+cunKfhQ69W5LteClH/w+fZ0X/k+gx+KzPyPATf2LasZpGMxzm5
fwmNQvCXn6W0LemB1PVi3r3keeL3LPrH5O2BjO3r++uayef3XLzvvLYHVntf
ltmR5zW4d35+fl0Ppi4pyP1soBBrDyrOLRXuge66TeoatRRizZ/wyY2SPaB+
jPpcSM6D2qS9DjHbehCVJ9ZVRs6L52dqmVLSPVBsKZYrLKQQfF+Xju6S68HG
6atqYQQ5vxzsB4+o9ODcXOjJPHIerRrU6ril14P04AuDlHgKUfmSdpZ9vAf/
5HW2niPnn5/mvVY3/R4E9RnFJUVSCJ7YtCZv4x78sPlov/U1hVh5fmNd1IUe
uAQXPXtNzlfuhvGiils9EFifnvvEiUKUPvLVNLjdg7n+QKVd5Hz2UpQqaHLp
wUiG+fpftuR8DzAjeh70YIVE0KZ5Swqx7Gjxt9lnPbjTqxFQcoZCLKG+SZWM
6IFnSo27jSqFaDugqdMd2QNDgeeLwsoUIvXbCC08pgdqZz5p/dhF+umrGj//
ux7cDpfe2y1NIT6ndDtxpPcgMOK0QBB5X5lHSKp2lpD7/wwLTGNyEUoilfVh
5T1wWEm7GcngIlaEOV03qupB55cC04f9XETGi6KIivoecBwRtBBv5yJ4fK0W
czp78MtXpLmuiLxPXZMLQ6d6MP9I/dZCGBcRMGt4/tRMD9a+8jyl8IKLuOL8
d3o1uwcv8mr9jPy5SP2Pb/dc6MGu+aTt9x9xERbXmN63VvTCe3a1/0FrLoLf
ZM9JfbFeeER2DMcqcxHXlb63rtDrxcE8vsiZUk6CqG38mH28FwHH6BndPzmJ
9bbj7lb6vfjFW3Mk7xsnQU3YJlVk3Av69Kl5s4/k+4HQc8eHF8n1QVsOigRy
Et/ZV3lmnXuxpySn9d0pToKHyqv2O7oXaeOJJ/bWcxDJRy9H5830gjG8EHzE
mnyfkpVQB7sX35qTGT4mi9Q8Xlo79U8vvvPM81K1Fqm1DeeFfnL0oa+dOcGW
WKSyzpsGFKzsw+MbRPD31gWquqP+/dLNfZDW/Dnqqr1ArQ1XP1Ov3QeDpRef
HV7/jzo9LLJmILQPA/rTDzb7zlMzDY9QQjR/48fc3Id1t1hUdx3PqIyB37jw
ke/YG98RqkSnbkucZz9OD+RztvLTqKoqHiHCuwew0U/DesnGr9Rzm2vYnA0D
6Kjr544uqAOvsGPm98YB5NvVapdJ1uPnuvUOzi0DUJt5o1LwpB5SK8yGhjoH
oOOrxiOt0QDm1EBTzdAALpSuc0mh/sLjsoUvbxcGMNCSR7mW2oz3TvI398kM
otpeYfNHq05Ml/v32j8axK0WyjFTwd/o4Gm89vfJIERFfBVzD/xG/rFNzGdP
B3F5wVaZcuU3nld/+BPpO4ie+hVH1FJ/Q6a+RKDq5SAuDMcPmGv3w7yV88iO
D4P4Xq9YrWE3gOJB5+Se+kEsO/ef0uyHIXyU/rHLtnEQlNedG9KrhvDCekn2
XPMgtrQuPXNifAgXh4ML+ToHEdCiGiKzl475sdROtcFBBHGNVdv9oENhhs4X
/mcQauObODeVMhC2xNzp1NYh1Bw2ot8uGsHCFaG9dVJDoJ9VSHvbN4KrBbXT
J3cMgUu+vCCecxSKjzXuHJcfwoSaWLCh2ihKF6Xv6vw3BP8EX6HdX0YxOz/l
rn5sCLF92zz6w8ZwesLbV8ZpCKdu3GylbRtHrr7GsaTbQ6hdvbaCsXscW9L+
rpJ2HULbw1zO3xjH2E27gG1uQ2iOrtoWbzoOzxGjIHGvIeQ12cec9RtHxqDY
K8GIIeRorDn3Zmwc67oyY7kKh3DgUUSYduIEAlY9eJBVPIT/THgkRNInsHK/
5pkbZUMgdjXvbf0xAa6w+tVN1UPo4ljBFG6cwKQB0yWpdQir9/HwJnBNorZI
+qQ+cwjuXl2ZnaaTOMZi7lg6OYSlyReTra5OonRL1pJvrCHYig7TOm5Oguqu
9U1yfgi9jmN/vJ5M4tO+y9tml9BR9UXZczB5En4pb/68FaHDIS2s5NvcJLg7
LjedEqVjS6fE0lnOKXislPm8XJyOhAulept5pnDPKtvSYRsdl1cn0uVEp2Aj
0VinqUjHscVRcYODU9B5yft+WIsODlOlwITHUygsaHwcrUOH2VUxnSCfKRyc
jDhvrEfGU6UNWQVPQeXkznVUfToO/fTpaIyZIv2u4xZsRofcm47rGdQpLH3w
yGifAx2v7rTcuDc7hYyJETb1Fh2qajekEv5N4arlmSjt23Q8Hd14O5fCQuFJ
BYbhPTr4D3Y1f+FjwWNLp/tNTzrGB9nvi6RYWFL6X3LcazpMB4xanA1YyDgQ
ry8TQcewW7GQ62kWrn7mm0mLpCO5V2Da2ozcL3zgYF4cHamf7/sLWLDwxDak
qSWFjlBu79P+d1jgWjfOyZtPx4sLjl68r1lIf2b2LqSQjlNh27VHIli4/Lf4
2MYSOm7Lrv6WFc1Cfv/bV9KVdORs/1kp9J6Fx9l6chpNdJjPVbf/zWRht1xW
fVkLHdHjL710cljojdniot9Ox+Zs+r9HuSwc9mXnn6fRsXVN1LG8fBY4zd+d
cWHQ8ft9SvG5ahbS6vkXFkbI/vn/ogjVsXDxiFucJ5OOJ7e5xKkNLBCKxmPB
LDr47Ic4W1tYsE8gXgjP0nFy6ZJq3XYWxEVk90WzyecTRj0/dLLwkELxSFmg
w3BrYvGeXhYUXW7uUOJk4I2T2enzv1mgjbRWf6MwIP0xudtpgIWDTWnCpdwM
pJ95O2rLYGFcdxP1BA8DlJZ3ridHWIimPrvauJqBT8fFeMXGWOBIMk/t5WdA
83zbBa8JFopHmm9GrWeg4VPBW40pFvwUDeTPbWBARvFwI53FglC2RkqjKAOD
C+u1Z2dZ6PiTeyNYnIHiOnkfczYLsVCWPSnJQFbjlc60eRbkS6STS6UZ8OKJ
KhH8xwJrZayN504GbA5utJFZYOHbiY0yh+UZEIyN2r1jkaw3OIS+oMjA6drN
29dxTEOraXXS9z0MfIiOOT5E8kqRZ9dclBmo2XDk43vOadSe55Deu48BVWvN
o0Zc0wiNcR2c2M/A9cHfMsMkn+ufTPykzsBebh8TW8o0tuy4YXn9EAMs0+TG
NpKHbvRvk9ZkINmtKXnPkml8SjPv/63NwL9tqbQ7JN9iNcfHHGVAajnz7juS
//vP4Kr5MQZa2ywe5ZO8cK9ccuNJBmaN2POlJBdSNfqaDUg9NyeMfSPZh/Ij
NsSIgRXVd6+Ekax/RPmygQkDHn0xFy+RLOibKsF7loErVJPR9SR3VEv3lJ8j
9XjLXPOVzC+WPzb62QUG1vjV1x8i+drpjRc1LzPw95q7ag5Zn3x4iBinBQMF
TT4nN5HM6lzd/cOKgVcqUcLXSX1yJJ5F3rVhYCDhfnQCqZ+bBYe5ii0DG0fn
hitIvTWTXEVZdgxYP16kdJH9WDk62ZHmyID7g6DxdrJftYo3ImydGeAJKsop
+stCqFO/mYwLA4cppVZv/rBwLtt84+BdBvQS27jOk/0egkH4xYf/1yNOIp70
x6cn5aaiTxiQuKfzdccMC04lGsJtngxcnG03fUP6afGEcpihLwOSX16NHyb9
VxicarImgIHXu89OuzBZ8GmSFqp6zkDMHci+GWVB0HxjqHYo6cfZ50bxdNJ/
MSHGlNcMyD8/dMZ3kPRf/+r1xBsGMsPfppv3syBnyxGyL4YBe/uy5+U0st77
/UGyHxmICk34l9VE1kuYG9A/MYDlDQTXL7JeSsuaxM8MnDD+vF2VPL9mvuWB
YlkMPGl7zm1XwYJjeKo/fz4D+Vtb90zmsZD1N/xZfiEDVMXwyvTvLPw1f/rI
sYT079OtPJeyWXgmed65vpLst/JYzNPPLER8XHk+uJmBxqwxo9lYFnp4Z04f
bmPgjupqVakoFqQcevQnOxi49ilESOMNOU/2Zmue6mVAvCtkQPsFC0V5FrL8
owwwsu0oLo/J/CUMpPKZDCjrWbUJubGg/+SAuOMkAy5KIpyJriy0HVknUD/L
wLKluTb+9iyM1RJ/griGkdO2f8z7PAsb+jZWrN0wjPgb7iMaSiyc11pe9FNk
GMEGqrLV8qSe7ybzHESHYXthk8uRHSzI3ihLr9syDMkLgxbLN7NwaPp2eJDc
MHTzHOd+LGPBZnm99VqNYQh8qNuxtnEKP3Z6ca+9OYwkrt2VHy2moPxLqCrY
fhh9rhun082nkHb/XZDArWFoMZzqk02mEFtVLLLBZRjN/C9bnY9OwctuqZzY
42FEvV9ieUN2Cqe+PDaQCx1G94YdlDHmJAZVH4QfzRvGRwc/4XDbSVzsW3Wh
ghgGp8KIB6/FJNp8IySPFwzjcaq20p1zk6js+P5Rv3QYCe3tL7bpTSLdbT7P
pH4Y8653K2a3T+J+/p0+i4FhVKmEDD7onsBaXUfZx7wjWLlaychBewLRP/l4
VNeO4H3Q+Qvp+yeguC+FPrFuBKrF0UmDihM4IT2UeEl4BE2vPLUUNk7AZ5m5
xKGtIziwWTxeiTmOJQVHBTlUR8BJD3KQfjmOmf0SHG4W/99vgss5nYkOudoG
19wR1Ow8GkH1G0FaJLb+Nh+FRbh9y1OZAZxPVV4rujAKLmWbOk7Hbjz0jHnU
8WEM+zvsLwTo12P1neW6l1LGUF5w5fnvLfUIv2bLP5A6hvmn9Xmm03X4orcv
jpkxhqXmOyRjw+swsLaukIs6hk5NF9v0gVroRXJw72gYA6uII0FoZQ0Esy4E
OP8Zw0eXz//1Z5ch7n3R6fl/Y3AaTIgxtCyDYvhOMXcOJnluPwWsESjD0Qez
qV5LmRDoCJy5Y1+K+xqBdW/4mEht0J5fKVeC3po8gXxJJjyPq7WnZRYiZUj0
De8xJiYEbFgW3AS83npe7DrBxGHx1l/e9lRcPTW67ZMBEyPlT7R21+VhU25u
2nETJn44jfQdefkDPoFmxX6XSd1K2Ns/SH/HNeXXEytdmfhDdRS59uArNIc5
struMbHTUrzYaiQTYtHX7n9wI9fTD6pbnstE04r/lut6MFGqZ/3oKjKg1dm0
yTuQida6+Jdea9OxxUNAZ3kiExptNxrnq5Ox8N/91c3vmeDYnd05U/UBraN9
9YnJTLgL9DEEa5Lw3OTLee3PTGwTlL9Q2fwOizKnbnnmMlGyY0z1OCWe/O4I
jKQ0kPmHMO9qe79G9tPZKw2NTDQ2CY71Vr1CyP4LO+JamIimHQjjEwqFXoJ8
xuEuJgJDNMdlsoLxzaWq7BGdCS4DScPC814IldsbaDDCxF0zjd8Tip5w6I0w
kmAy0ZNlZ+Sx+jGkj93oJlhkPYtRvpNFd7CE81f881kmVM+YV7VkOYCWud/m
4jwTwi7Wr0JdrZBrE6eg+I8JMxO7ImtDE7wS45leXGSCtvtTauaL2+r/A13Q
BNY=
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{
    FormBox["x", TraditionalForm], 
    FormBox[
     RowBox[{"u", "(", "x", ")"}], TraditionalForm]},
  AxesOrigin->{0, 0},
  PlotRange->{{0, 1}, {-2.494931344313982, 2.4970526788478478`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{3.504351096388619*^9}],

Cell[BoxData["\<\"conv_diff_u.pdf\"\>"], "Output",
 CellChangeTimes->{3.504351097461491*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"source", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"u", "[", 
     RowBox[{"x", ",", "a"}], "]"}], "*", 
    RowBox[{"D", "[", 
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"x", ",", "a"}], "]"}], ",", "x"}], "]"}]}], "-", 
   RowBox[{"D", "[", 
    RowBox[{
     RowBox[{"u", "[", 
      RowBox[{"x", ",", "a"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "2"}], "}"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5021991167832737`*^9, 3.502199165625288*^9}, {
   3.502199445286195*^9, 3.502199449171276*^9}, {3.50219980583748*^9, 
   3.502199871761393*^9}, 3.502199941959586*^9, {3.50220002712248*^9, 
   3.5022000282904387`*^9}, 3.5022000861529913`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"20", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"a", " ", "x", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}], "+", 
     RowBox[{"Sin", "[", 
      RowBox[{"a", " ", "x"}], "]"}]}], ")"}]}], "+", 
  RowBox[{"10", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "x"}], ")"}], " ", "x", " ", 
   RowBox[{"Sin", "[", 
    RowBox[{"a", " ", "x"}], "]"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"10", " ", "a", " ", 
      RowBox[{"(", 
       RowBox[{"1", "-", "x"}], ")"}], " ", "x", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}], "+", 
     RowBox[{"10", " ", 
      RowBox[{"(", 
       RowBox[{"1", "-", "x"}], ")"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}], "-", 
     RowBox[{"10", " ", "x", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}]}], ")"}]}], "-", 
  RowBox[{"10", " ", 
   RowBox[{"(", 
    RowBox[{"1", "-", "x"}], ")"}], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"2", " ", "a", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}], "-", 
     RowBox[{
      SuperscriptBox["a", "2"], " ", "x", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"a", " ", "x"}], "]"}]}]}], ")"}]}]}]], "Output",
 CellChangeTimes->{3.50435109765758*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"J", "[", "a_", "]"}], "=", 
   RowBox[{"Integrate", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"u", "[", 
       RowBox[{"x", ",", "a"}], "]"}], "^", "2"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "1"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"J", "[", "a", "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"a", ",", "amin", ",", "amax"}], "}"}], ",", 
   RowBox[{"PlotStyle", "\[Rule]", "Thick"}], ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"a", ",", "J"}], "}"}]}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"Export", "[", 
  RowBox[{"\"\<conv_diff_J.pdf\>\"", ",", "%"}], "]"}]}], "Input",
 CellChangeTimes->{{3.5021991167832737`*^9, 3.502199165625288*^9}, {
   3.502199445286195*^9, 3.502199449171276*^9}, {3.50219980583748*^9, 
   3.502199871761393*^9}, 3.502199941959586*^9, {3.502200159339991*^9, 
   3.502200160153182*^9}, {3.502200345115892*^9, 3.502200346856551*^9}, {
   3.502200629524897*^9, 3.502200629849448*^9}, {3.502200887304186*^9, 
   3.5022008924055357`*^9}, {3.502240917621992*^9, 3.502240921581522*^9}, {
   3.502241047689166*^9, 3.502241048431182*^9}, {3.5040002760144978`*^9, 
   3.504000286004998*^9}, {3.50400086781607*^9, 3.504000877470779*^9}, {
   3.504001083155376*^9, 3.504001112458433*^9}, 3.504001198935052*^9, {
   3.504001430705776*^9, 3.5040014649677343`*^9}, {3.504001497503985*^9, 
   3.504001511549205*^9}, {3.504001554989911*^9, 3.504001596705283*^9}, {
   3.504001676933679*^9, 3.504001683474367*^9}, {3.5040017315958967`*^9, 
   3.50400173223351*^9}, {3.504002308508747*^9, 3.504002325570188*^9}, {
   3.5043509211648283`*^9, 3.5043509794399853`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {Hue[0.67, 0.6, 0.6], Thickness[Large], LineBox[CompressedData["
1:eJwV2Xk8VF8UAHDM2HfGNsQwKCHJkhT3SAsKpbRIhbJmSSRLJGTLVsgSyq4S
SgotWmQtiiL7LiTbT8kb8bv9w+f74b377nnnnnPmM7J2Hhb2LExMTJP4x7/f
3+ZN1ZmY5pDwnArfVNgUquMrjTXjnkOiPeWb1fOnULmywFSG6Bzyc402vvp+
CsXZt+VsU51DayMcxHa2H8iw6wjFy2oOzXTlp4pH/kAlr08vjT2ZQ/+13Gs/
FjuNrsV7vmx2nkfncqlDuzJn0BbVZOOU9gWU8knSKNthHsna2j/iHv2FFLOl
F4wYv9DzeLJQ2Os/yOWere82OQLl8bJpZMUxkIdJhxLx/S8K2tog0JTOQPZ8
5+x9lv6iY7bRM78KGEhXTkJ1im0V8VTw3TetYSCRK4FalQqryNtaTPbvLANt
zxVcnbZbRbsebOC3PriCXgeMqpf1rKIRo30/JCh/ETIRLLGsW0NyYQm5yWmr
qHtp4gW5hgm+T2rdmM9bRdOtDhTp90xQbN4TZFq2itJ0yqu2NDOBlqSiFWv9
KtLYrDFk2skERo+fC1xcXEWOj/eZ7ptlArfB8cDD5mtI8/uF36elmcE6S7o4
ioUJ6mzyg3p9mUH/4HzOVg4myGM+OHg8iBlo5Nq0MV4mWC1x02oLZYYxZ+cI
Awm8ritXWXkcM7hrPbFbVmOCtBDnnG15zHD1g7GE80kmGCxtdDZvYYYCwjvc
qJIJ/u47xZkozQL+NeNbz79igmOsFQuf6SxgGnZsMqWWCTxdTDO4lFhgkXfH
/u+fmIAjnGW9swYL7JQlCYVPMoHD2yM93XtZoH/vjcx3VGZI4Crtj/NgAdHk
h0/0A5lhqX/sg3I1C0wdl3FwwM+ZotLWkvWKBV5K3xCLi2IGtSaKJfc7Fjhb
5O3fl8wMjfb8k5+aWeDx8+1wuYQZFhipnHx9LGA23NhcNcAMfU9VJM1XWeCa
2viw5k4WCJLZ7T22nQR/Kr91IiMW8KrLWrqISOCys/mDiRkLeI4xSTMZksDM
suypjRULTHOKWDCbkEA0wD865gILPAjx2v/sKAkK6vi2jOawQDlz3d55TxKI
H2BeP1vEAqMi/7VLXiRBdNd/kkQJC7z9xbOGfEngPv2NVeA5Cxi8mVS/GESC
rUK5XdvbWeAz70+uxGgSNJzcejWRRIK9EdN5Gdkk0P2u5JPFSQK3JdmNtnkk
KD4vde4ePwme+jcfkCkkQUIos2WNJAkO5jGHRRaT4Pi95g0/NEigPisSx/eM
BFOLpz/tPIuvCxktmGwigXWQxXtTFxJ89IkKMvpIghb23dXHzpNguegeU3Yr
CR5TN+a5XyaBqVIgO/pCggBYvJSeRII37c+11fpIwBMTKbNQS4LkwKipwZ8k
8C53lUrG9w2WKGGrmiVBT/cBCZ1PJIjS0nt4fZ4E95UkhIN6SCD7xiGT9osE
RvX32LgWSDBjHajAuUKCUJbmn7IyZBCtCHv9lYMMU0qlU7XyZOg5WLN6hIsM
FgcTvztuJMP22aJvbdxkoN09MfRQiwxaJ/3OVvOR4ZXe9Jdt+8kQyXd7zYhC
BsKX98VBPzJ83NXhbi9NBru785WLQWQgvTnJegOv01j/tSIljAxLWUdNn9HI
kCqaVdoXT4aKjOKl/+TIoP1kU65zARnIbsJDquvJkNUjdJenmAxyjns3qW0g
AytpKaP0ERm4Yzy3KSuR4cvBmlu/XpCBS2V7grAyGS7MmkdfbSdDrX+ie/4m
MnSLakbId5HhUJ1J92U1Mhjoi4fV95NhnP+6oPlmMgjEDgbyTpHhkVrN5Ig6
GXyf1PqXzeL9ZhnHFWwhw0BP0aVDv7C1olnPapChZKOnZxoTK+QpRHq1auK4
WVi672BjBe32FbdAHJdAv23nBrhZQX/frZ2K2mTY38BsryDGChyf2q46bSXD
pP3NY49UWeFY8o240m1k2HlHRKBFgxU+fl2JQbpkSP+WWj+1jRXCtpxNaMTe
v++OtvweVthx7kpZ83Yy5IXJzcB+ViifcG0z3EGGvy/z809asIL180bSU+yS
zcWUlFOs8Kvtyp1IPTKwuah9KD/LCuse6PFPYp/OfRz6yYUVShc10gz1yVDZ
q7X953lW0NwhB6nYAqJVC5yXWGHI4iPbd2wn8x33FQNZIdh4bVYNkeFNZI2t
Yei/54lYvYBNfbtTwiaKFZx1dLXKsC8w3n+6HM8K5rbzt75jN2kaR6Yls8Ls
6+uKEkAGuvsH9PQ2K5Qd/Tm+C/tyoflSWzYrqO5Z63bB/jLYVjJbyAqd9aWs
17FVqUcceEpYIV91wTkfO/xQ1zqlJ6zQ9O0VVzV2f4z1193VrGCgJfyjAXtr
3UCM3WtWuPZknrcNO2HNbteVOrz/Blufr9gTOuOM2x9YQavXY307tsEF5/LK
Nry+hap8M3b6g2mXr99YYUNAmcdL7IVRD7mFflZIvEnhvY+9T/q/Lr4xViio
9FtLwM476nND+QcrjE2TTS5gryQsGxnNs0IKx+CsKfbhpstM9ks4XlXbFunY
D0lMlVf/ssKK/26rRRwfNr1QjywSG4glbleowT7lw7b+OScbmHEdPxaG/aw0
qr+Tnw12FA78NsQWmOS5tSjCBq8mRMir+P28OSHMukmODY59eHvJFpuafOuF
yQY2CG/rG+XCvtAi4e24iQ0kdOc/luD3L29AG7mrywZDrZdlx3G+XPbPTX8J
bDABr276YX8tV7To3sMGF2KjY9mxwxVV3wgfYgNbqUhVYZx//afLfDcfZwPd
hoKuBJyfW9M0NpuexteLKq/jxJ7g1r0Tfo4NJk3P+07q4Pjufnkk15MNdvHs
TDmInR4EfK8vsUEd7+0TT/B52De/J3A5lA1c7TRHnPF5ydvYpCkazQbyO4oe
vsDnaeWM6fSWBBwPKy4VLuySjkMnXDPY4JqR2qFkfB4FXthsG3jCBgPqTI9d
8Xk+lfj1T2I1G7T77/+Ti897sYtJpdFrNtircai6A9cDIwnNrY+b8f7/63y5
SZUMwT7sWteG2YDfxVbiHq4nLaaXF7dNsEGiJl/SC1xvpBTmy2d+ssE51cO5
zbgeVbZ3bzm2zAY2FYdiehXIMLu5ZLOyADtkXrYeapQlww4O+tyACDv0H18+
Xo3rW/RASmmSJDssFzcFF+D6pxh3ddOqIjv8deYb91yHn/fHYZXPeuxwI6rV
5as4Xj+fsf7SOXZgGfxo2Ibrq1Tg+e/KnuxgzWhQvsJLBpfDYwWDPuywwd5g
aT0Pzh9Sq4JJCDvEfj66x42TDPo2OfR1aexAOuugUE/C8aIay7yrYweRHWfm
q5dIEBt/S5RfjgP8U0tjn/eSwF85RqJwPQe87/XtWewmgUP9VSmkygGd+VK5
G7tIAGuuch46HMCza+Th9a8kWHA33NRqxgFzyFpYvoUER8zmdscHcEDZsqTG
o1ckoPGaXOT/ygGjxtfeWWbifnUP+RZ2c8DkWW4Vhdu4z+3SCkCDHFB08W/w
fCoJ2i7Trnr84AAH+W+rwbjPhU3/jmll5gTqk/hHAddJMNmcmxe/iRPyvitt
HsH9+3H033b+SE74rPSlPOkg7r86N2/IxnHCwaEVy89muL+PK5prJHHCkbyt
41z7SaC080DTkbuc8HVt36j3HhLcXc55nVnJCSF7SAPyeP6IczJ6qDzJCbNi
W4SW6Xj+2JMYvteECwI3X3yc9B8LyLIobbvCww09zNEni6+zwPjv5+lzPjzw
6fRyVeUzZvjlKJz3aY0Xhs2MG3Pr8BzHlBF79So/LBTph9e+WkPW3HXb0r0F
oDH5Id1xyyq66j10ueCSAPQcctgvprKKCvpWah77C4D/8UqLZjzPzpZq7GkO
FgBa4QMuJLGKQg5nW6zECEDq65RCS6ZVVJh5+dypfAHw/B7va9j6Fy2oaWTK
dQjATPpCl/75vyjy8F2mB1sFgfXqRZ2wmhVUFPNkb6iuIJS+rp7JrV5B9bUN
cSf0BKEhRPVDfcUKYtOal+IxxP9PHiBLFa+gUJGd29zMBCGaaSL/d+oKCuoY
8VSzFwRH9uuWPl4ryPvY+pHym4KwI+hL6gPlFZSUsH1jTLIgfKsPeHVo/Qoq
bzD3PJsqCLbj5w2Y6CtoQcd3jZIlCGJbHM86UFeQh0SjpM89QchRXGfgyLmC
XLpdDuu8FgQyz9EjOyYY6LR1yfsXPwWBu6GuQvceA/msE5r6PCcI6DGUyOQz
UFz/Rb7v/wnCoI3BJGc2A72y0TsqSAiCQUuo0kwqA0nbf5iwZxOCxz5Mo78i
GajfbYqbX0YIOgY4SlhcGOjkFUUL2wNCoHTG+ZTwZgbyhmgfn0NCwB91wDpb
hYFimGfSrx8RgsVORo6mEgO9CK0YqbAWgnyrU5zn5BhIKnLXRS5nIUhtmVAT
E2Gg3gS71PKrQhAsre/ryiDQiZysftZyIWANbxVh+0Cgt+/f5vE8FQKviMhl
o0YCbZwcdxGuEoKSPZT+hDoC/VFTW6LVCMGlGK58tTcESn5Zw7+jWQg+Xevl
uPeUQK3fBuHCiBDsP9PCFp1DIO0VMrvfuBBsUym1kbpLoCwZpY9XJoWgaati
QXkmgVwdPI/HzgoBZaLs51wqgTgXmS8UMYRggq9/Z0E8gQz56bn9wsLgFH8i
nggk0P0te53HRIXBNW+grzuAQEJHzqlNSwjDeWrFSo0fgQYzyp8vywiD5Upl
bupFAgVu3PWFoiIM4yHFNqFuBHq225513y5hEGe2/bTjFIGULxc6PvMWBsMy
Mf3NhgSKca9Vj74kDEIPY+tfGxBoxmaIYe0vDHzvGg0tgUCPd0smsAQLg3q6
MU+0HoF0+eMrTa8Lw730mzlyOgQyyfHhHM0WhuiEn8x8mwhUnJTY/jRPGFpH
M5/1qhCIL6IsM6pQGCK+mOqVKhOo7dyUutpDYWBSDSk9rUSg41qnrPwrhSEq
ymrmpzyBnBp2PxBoFYbILs7al1IEaqq2uzjyWRhyDvw53SBJIJWHV9DTL8JA
/lqb30El0NyNqvYT3cJwqoVZYVWcQJdOqK4UjAnDbHvhA18RAn0zNanzmxAG
fR+dA3cp+PnBMWH/D2HYvd0tulmYQKvy2fILc8KwkuFeoSJEoPAZiqneijD8
FzOUIcFPoIlBdTGBNWHwdajMP8mH99tuNjTMTAGzaVDL58X7q4y8GMlOAWb1
qip9HgLdurqS2S5EAesPsza3OXH+eEk4FYhQ4IyLgDyDg0BWDtpb/MQpwPnb
wOIUtvS+83XS0hRY+PV1dBM7gfIpozPOShRo2L0++Q+ZQEwJB275qVBg26s5
CydsK+6XelFq2B+niB4SgfhZbsUUaVHgCcHCaGEhkEsgi2alDgUaBYf3mWHX
/nHvqd9OAb3XXr6fmQnkN7t343cDCgSH8B4ZYSLQF5cnn3/vwtezvVzzxN40
TvNjM6JAiZC7Dxl7pHe5XsGMAleKmF9ori0j06YH9metKNCRHYze/l1GRbvF
eb1PUuDtz5EJd2yWN2FPQm0o0G5gaSCD/ezZSVKuAwXMLdXnoleWkZBG0/3H
zhR41LiivwfbtUTb4q0rBVpqUlRZsWXz+bOHLlBg7d5AdyxjGV2mXTaav0iB
/U86bxzF7rg9Mcvkh+N5WatLHnuzmGWKwGUKuHKJP/xFLKPrN9/o065Q4FvX
bc4m7HHeTeNqIRQwHOqeysaGqPRYdI0CuunT+wOx08nsWuaRFKByjitbY/+6
4tV76joFumu7QvWxDzAGQt3jKDC22nlAHvu+z37loBsUmI6eiePFJi9UtsUm
4f2Mb9AllpfRaTcF/8wUCvQZJJhNYVdP3JB9mE6BL9eUG/qwKWdXG15kUsDv
5FrOF2z3AZfzH+5SgFePr6cFu9GqU6w3F+eD8zHXD9jyHYY1PwooEMDSa/ER
O+hgmQPjHgWqWJOjP2N/+yDFx/2QAvtOBPJ3YWsYRVVQyygg1xPTN4Id9+6X
9cZyChhbVs3NY0/o25F1n1LgTdrvXST8/DurWx4YV1FAJVi7Qww7U2v7oeMv
KFAw63RXDftPWSHhVEOByPGAOybYFiqUHN+3FNgYcqLVCbu4MNg48j0FctkW
lKOx2ek/51IaKHBzfP2jEmy7rOOphc14/00jVh3YLyTq0LMW/D7TuRSY8fvz
FLgT1/GFAj0EncsG+8N1bu3xTgrYfu2USsJWZPft+9VNAWc9P9SM3fP3gIrI
EAWm7H/l7sT5ou33sl1+lALjSWF9IdgJi0oBmt8p0KvnIlqHvfsHS9OhnxS4
Y3LjoiXOx3Qtp4tqcxRYDzc9MrDnrnykcf9Hga0San/G/v2dknrp7R8KULwy
CoNXl9GMnqqCOlkEZNseynzB+W8YefMTD7sItJ/8u0kNn4/UtqWACU4RyPCN
fxSDbej4tu0OvwiU7nnJYo7PW2rC0WA+SRFY3dXvPYfP48/uF8pT60SAiPj+
8hQ+vzsV5Drf00Tg9s6cU63Y01XTqoGKIkA/bpn2DJ93g5GrPT/URWBv/MOb
uWwEmtIq1m40EgGFl7zXx7kIhIIFh/P2iUCJSE2JLzeBkpp8YoPNRMByU+gC
D65H6LTB6NbDIuB4Zt+H7bheJUV23Cg4LQJm19YcHuP6ptfD9DPURwRE/6gr
MeF6mRBsmaufJwLXI/9ascoSSFM1YI67UAQ0LPfmfMLu6rqr13VPBNZO90Zm
yhFITmO606tUBHTzqhR34Pr9ZCyE5161COgd4DKOWY/rr0nZReHPIvCdYbv3
Iu4PMiLcRhN/RaDg95Udh3A/Kbn3evrGUVHQC864H4X7074BLbEMK1HQqjjK
//w0rr+UBwYFJ0VB+MFn5582eL3g5FvVZ0QhvmLj3OEzuN4ecYERD1Goeb7Q
vsmJQFdIlCTNSFHIVlqI4rpAoAMnHbd3VorCs2OfC/eEE2hBgC9aiioG50/1
FTOVEWhIO16hqFUMBlTyTmj9JdCF2v3qz9vE4JpE6g2fVQKxWHDptXwVA0de
GYlna3h9t2uHF3vE4IGqrJA2CwOdzQ0IhQkxqHnkaKvMzkDj/E5D35jEgc+V
RXVegIF+TEAm5xZxeLexkOUbnYEu+6wWSWmJgwljLGBVnoF4yC+eqOmIQ5Xd
2BJdkYFUZLZ+sNQXhzmqxey5DQzkbqnKyDYRBze3paxZPP/MvZE4rntGHLym
NdbXajHQr/QFyrkkcZBjb+pN3MNAX6LPv/VIEYfcB+4rN/Yy0BP/GQ/vdHHo
yD/6I84Iz1NWU02Bd8XBuCVs8zUTBlqQGA5OKBYHGVuSyhkzBppN/TxdUSsO
d+48vzFiyUCTt8pqmX6Jg6OB9ebpMwzUEL7pAusfcXA4+tfmzVkGKvIpluFi
iENjS4xmsj0DOR0p8hdmloDM4MkaHUccD5G7mxV5JcBEMN7IA893I0nxGfsU
JOC7aIRW2nkG6rvp7n3LUgLKlYyGuwIYqPLgmILjMQmo03i71+cyA90UsO7c
ekIC/MykzgoGMtCeOGPdLhsJePwU1RgGMVBZlDyTlKsERPSFz6UGM1BIcPf1
7BAJCE7WbiNfYyBrdFDvwjUJqDbPd47B1l6tn9kZKQG2V9zLhMLx+wmoODga
i9ffEGciGcFAhy8liK1PkwCJ71yFUlEMpOi+J7e4VAK0B9iqB2IYiEn15eHA
xxKwU99VxiKWgbp/aLCZVUhA00YrtXfY8c6yzrPVEpC4uGSaHcdAf86ubNpS
JwFWz78HmyYwUBvda5ClUQLao5W1KrCLhydvtDdLwA+zmJuSNxjI1qZz0fuz
BJA0aUaj2E0nHlc/65XA/eXqbftEBsqjKrlGDkhA2qPryW+xg7rurDs+LAEq
UtNbpZMYSONobDDxXQJqhYIufsbOsHDao/efBDAXC7Svv8VAHsbSbW85qeAi
4CeRhOfxnTdGbT7xUOFrjNLNIWxK1/3ZPn4qHNpKvaSaxkDVTlt5l0Wo8PrZ
3vDX2GyRB4zU5KigIWvzqj2dgbpaRTt3KFDBYn1EJuU23o9Yn73JBiqw7bWe
OoRtUegcar+JCiVJbC0fsTPrQl7d1qWCuHlTbUkGA53nMza9p0cFHaJKfAR7
5xH+3qdAhbtPU3pFMhloYuz28uc9VAiL8K/wwX6uahc5YEyFNpp2cx523MUN
Yj/3U2FP636Lz9iarBWaHIeoYKsm+VQhi4HYTQPeiRyhQkBqTqIpdneSgQX9
OH6eK/KjXtjB8i0e+qepUOHRNfAc+5Br0uo+OyrMk35d78NWfGIVe9yeCqLr
e4v/Yjcbfr/vfY4K06NrZtvuMFDW9YfbQtyp4BCk9/kw9vl2r4Z4TypoZRk1
uWMbSuoezfSmgq75Pd0IbNEzTOP3L1Hh+O0puSzsyft13pX+VLCMTA8tx36x
EEOqC6RCMcXkZD12vO6hm+3BVHDTjC3rwrYLkZAdCqWCYAsROIWt1TRQOhNO
BeNnKnXL2BxCBforUVR4+6gmiv0uA/Ucd/3IGUsFZqdDjcLYJdlbrMUS8H5c
oq9JY1+d/DMln0iFcxu5atZjH1av8dtyiwrBhL+PGrai3zUOSKNCNntksRb2
8ut9KaYZVLh/dMBWF/sDh5DiiTtUKHuvnqKHnXXg2xOnHCpUrm4xRdieqVmG
PvlUcA1KC/vnXYNn20KLqFB1XVpfH1tsg7LtjQdUCNIK8t+OPeUxP5tV8i9e
l3S3Yr989iyo+BEV0t2ag9SxE9YCeaufUGHDiL7RRuwze3dl1D+jwvCL1BRZ
bO14LuWv1Tg/XuW4iWFzdn6qGn5Jhds2Ss3c2L3SKUZzr6kg18QoW8XxKnU4
2fn3HRXIn37IzGGHlNAduOupoDn3UXwQe71+WajiRyqwP6p/8Bx78bCHKNOn
f/E/UFKE/fbcpntdbVQ46vnVLQn7VFpxS8w3KqyeL1N0wlYpO2fr0EOFtF0B
sgewl+s2LqJ+KvRGbOzVxk5aLJL4bwTHO10zienf++V2Kv4wToWDh4QTR3C+
qcmtRwWT+DwY9x9+j91knn/2+CzOr2xL+WvYTA+yS2sIKrDeY3D8xvnf8sZm
Z9pfKtg/sNFqxb79TebrBSZJmL0vqVyIrc2WRSiwScL4oR8XD2G72qbvjhGU
hKGLz+9m4vOn63v8mz1FEnY1xaq5YbPHi59DYpKg9yYxbjt2zotbCQtSknC7
rOn1F3yev4kl9hzbIAmeS5r8y7geFGyycN+iLAnuyRM6NdheuwWZeTZJgrmi
lG4YNp9XvGKNhiRc+dtQx4W9q+W6pwKShKdxf+lcuL48CgtjX7CUhCMxyo1t
uF4F3TZMbz4mCWrMpwbCsfc/ZlHNP4Gft3ygdxv2xECwxTFbSZiQJW5kJDOQ
9PbAzFeuksDf8JfZEte/6PmLGtdDJSF5LGjyJq6nhS6/HI+ES0LoUn7zVuza
Ee8M2ShJyJvML+vB9Xj1qxdrZZwk+Gl6RNCwPas9O0fTJKGcJhv+73vQY6Fu
AfplkhDglvrc9ToDKYjYv5vvxWb72Hs8DNevuNGlFwOS8O1Sjt5QKM4H9rMq
kcOSUBlk/d4RO3XZLll6QhJuDhuYXQjB/b/fxnHff5KwOe6EYyDuXwsF1tz5
nFKg/rTtvSvuf691Dh86pi0Fi0zxXyq9GKhDl25vs00KHAOoHJLYP3Ys+Djt
kAK/5Dculy/g+mKQkO67Uwr4Hsy82+GJ36/xh6FUMylgZL/9r9wdnzerXee7
HKSghjn5eaAzni8CNGOtUvD1h1FW9CkGknpFqbf+IwUv8vNLJvA8QWfh3Gb3
cB0UHeNd2CzGQOZhA6tZntIwdF/d0WmEQBYPZ38FeEtDx7rBWyXDBDrSsTZ9
7JI0pE0tSPwaIpD1BlqPUKA0kOKzpUMH8efTDzbPwiOlwSbSHnL7CBRGGfZw
vyMNVYuV7YudBKrOGx3U+ygNUkWa50Y/EOjVx8VO6idpGGfdz7sV+81vcutS
mzRMTx6hRzcTqMFI4eWjb9Kw+zS5VL2JQF+n7dMURqXBo/fOqfB6As1pfj/I
y5AG/7ioyT1vCaT4fvJt7wYZGORJeLGxkkC+B6pPPFCWgbHSuwUJzwjU3BO9
6LdJBvS67KlLTwl0fkFZUUxTBt45zGS/ryDQcxn3qINIBkL2t+9yLsfzpf+C
eZ2lDFRZnvNuLiFQgBqjtzREBua/ZHMwFRCopbr5YtA1GRBiX/T2zScQbU8G
n2mkDCRYeF+ZzyNQnbUeTMXKwM1CDc+xXAIJRF/JlU+TgVcJrR1fsglUMEpy
SSuVAf3N3DJtmQT6lMazFNIrA5kD19eMbxFIskTFLW5ABmxDp30+JhPI4e3+
kbRh/PfnktUHsf9OxbSWTchA1v72fOskAm3cwVvUvygDsauVMpdv4s8fs+Eq
mWw0MBi+uzYURyCn5Vvf5Tho4DTyxvkCdiypMKeIkwbSGyMLSNjfxOrFn/DQ
YID/UtqGWAJ5ADtrsxANONqT3C9fJ1DGjYi+P9I0CDc2WLc/kkBvb6ekBtFo
EL1nU9dkBIG+5xceIsvRgN/ioWsk9pbq+kZ+BRpMG1pJNuD5vHGY/amiMg24
brdlHrxGoN8akXGHt9KgKWjD2dgQAknppxp369DgcJn3Yx1sA6Miso0uDVpj
T/SMXiVQjHWD/zk9GqgK+79E2PRrHI4hhjToEhOvW7uC31dHJCozp0HcxheM
okAC/bc3eUXyIA1Ss2lPHLBvVWVXRVjQQLPyt7ICdm9GtcZpSxo8W2C2yruM
43Nmej3fCRrwuseXPgggUPCcuYCrAw3aDSe39/vh9eysP3Y60iBruWGyELuu
3Sna0JkGmV35Ry5g8zy7SpZ0pUGyELcbJ3ZqYPmfRk8auFkoRIMvgUq5xIbW
B9KA9P6Ee5sPPl+X6VmJQTSobXZjKsJe/Kl2Yu0KDb5vqDAPwtb9bPS1IwTv
9/0uBVXsuhT/xmuRNEjIrStMvEggZ86I8PkoGrx0+yHtgc0TkGh48jpez2v8
4D5si1PFLzXjaODYljLFit0v3/9oJJEGzjObhkK9CXT11pS7eTKOzxlLdXts
eY4l5ee3aHDLIHf7XmyXH/wFN9No8COuLpAX+/cjg3SDOzQw3mL7KMeLQGl0
s6MP79JAV5f5QxT2jmQrikQODbiFP2R5Yof4esXN5eH77f54zhCbD+WH3LlP
A6W9927P4s9vqhePlygV43zU9unrxt7/gLe7/CENvqD2xTrs62I+6g1lNFD/
U5yTjX3fdONJi8c00LA7qZOA3RTaH9lbToO+06w5V7A55vYMzj2lgcWJm7M2
2OsVGdwBlTQQy8hqs8DeY126lbUa58vyXMRubIebZ87EP6eBsr4E/zbsaw1i
8RIvaSBbnO+oip232lyd+4oGVBb6DTnsWs3gcdXXNPAOkLoujj3ioilU+YYG
lm2SlvzYLNkTejvf0cDT4N4MG7ZsZ4bzh1oavL+27+iaJ4GA92DykToanBuM
iPuDbWPI+mawngYnxYcTF7Cv+FVNuzTSYO/lXoef2FmlbuK/mmhwaHSJbRL7
1ZjsrisfaLBSXu41ht0n2eHB2UKD3p7Y4mHslYPRtxNbaRDbw/JoEFsyUr9+
3WcanH8XGzyArftqYaGwDeendhH1n60WC6S3fMF2bQr6Z7+NJ0xefKWBrca5
4n/Xp9rw++zpxO93fi3n3/0rb73L/vSNBgEFK7b/1u/8cOmjVTc+b7PtExPY
v1lUlkd78HllW9L+9/wi2wblPfpocLru675/+9P0SDqw3E8D7fv9iv/2fzjf
6HLoIA2mPgfVr2J79awU8g7T4AGLnOq/+N0UfNSeMkKDLfe2HOXDfrzXfk12
jAaPeyR2i2HPlX88oj1BgxupnC7K2PxTV0NeT9LApKzvjjb2Jpp2ickPGty7
IJOyE9stJov19AwNoiq9662xY95aqE/O0qDljPGaC3bxH7aTXvM0AF6+337Y
U2c9nkQs0mB3Uc+6dGzO2/RBwd806Lh92/QB9obPndwZSzTISa7a9hLbUQ/O
lBE0qPrtsWsEO9xrMU53hQZj/hW2S9gF94qqa//ieBqW6/Dg8zAqIijUxSQL
r6LlBXWwSfvf651hkQWt0O2C5thyIX7OP0myEPT0YIMDtu3M0GtmdlkobK44
lordX/fYYyOvLJgPP5oh4fOpsGmv60M+WeCpHuOmYbsm9zhtFpCFA1e5vuzA
Js6Q7bSFZUGn4aydL7YI85HDOyVkIa/vy5klbGunqQPvqbJgLS7wXhzXi5zW
IFMjKVnY+iezWRd7c1bhHjMZWVD61Po+GHv/9mUdKwVZWC5RsqLg+pSYE6vV
pygLTk+yPutid3PKbbHZIAvlytETtthO30yUHZRlofNT7fwj7LCLGesuqMvC
/sY06+OXcP8tRSzRerLwu//C91pcL1nEvqzyIPw8KzHEIrZRkBMjHmShwer5
XQVcXzv231hMNpQF4ycyuZHYC5PD49nGsiAl4uRs6Y/7oXx4c5WlLMTsuLWL
E9f31NQPSZNuspAvtu/IF9wvrB1PJlp4yELVau9ngWDcv7Vnbjw/Lwsf7NYq
TLHvtfPHx3rheDmIVzRgV/MdilL3kwVnUz//d7j/9IZ1BfqFyYLbxSdXmkPx
9Z7jDpy38f0K3x6g4f5YZMy0TakB74feRi2/QaCJXp2CskZZ4FsWq+LC/Vnp
vKewTrMshPsPpdlh3781/HNPC46fpYiMcCLOx5HanLNfZIGzPkXqEu7vZYFR
PHcHZWH09DsOkxQCVT0SGhRdxuu77mqh4PmhSUIhgrxRDmaqP7Sy38P9v/HM
YZtrciD6UGZCvwbn6y/fpLpwOSB367bfxOaVjfuiEikHh2avpYxjf/OtPLwc
LQdp/Sb3Yl8TyH0Dr+XNBDlQUuAc6npDoPSIp5a1t+XgXDmXjEst7q+7OY9u
eCwHoXm7hu0b8fWe0ilx5XJgWXOzsxz7VaZG5+ITOQh7qJrEjOe36F8nj755
JgdzXVttM7Bp+Y+PWr2UA+9zhjta8by3n2x9LKZBDqLj068pteD6+7bk+Fy/
HFh3Hic/bCMQW/k5jouDciDc7tL7B9spd8Oz5SE5MLR4fnpXO4GUw3IopDE5
sFWasO9u//f91q1WkR9yEKi415blK54vGwJ3b1+Sg6gz8lkGeB7tb9mvHs5P
h4WAdc/8enB9ruEc5BKkw9JJMlGFnVNaFxcvRIfLpXHflrHtE2A6VYQO/BxD
Dr69BJo+qFn4QJIO5n1ML8/jeZfxRXLd5/V0kE7r/206QCDxnikOKaBD9uSC
WCWen5e25FWFG9DhLM9gyXfsjuiTLvM76UA/kDQqiuftJN1PzXW76ZCY7Kjs
jS2YXhHnuY8OebmBtI2jBOI6HkypP0IHSWczv9AxAk2WbXuvfowOp1p/DBdj
N3D8dzHjOB30TNgnvmJHVNp3elrTIfY+qVZxnEBk8X3pUnZ0qLMzf/UO+2+H
CO2CGx14HvWnTn3H+byp9VOvOx0OjEwW8E4Q6EV45NW95+nA2M1sthnbX5sx
LOWF7+/71cEb+3fyYH69Lx10jjTULGHPHX6gvC6MDjYiuktTkwRqLT7bG3GN
Dq9HeyLZpwhUQpaOXQing+/Pn4/o2K5PEmbqo+iw77F4kxX2pLDPowvxdLh5
6pFALfZwG2xtSKeD1ZmRqZAfuB9u2Cd6MoMOcYI/vW9jHw+y/DWfSYfAytKA
cuxPSi7lktl0CC7fuTaE/So4Uc2jkA5D7eL9OtP4+Tsy+cj36DD+08rUFFtL
pehn6n06rGPdsc0Wu7jzxYO3D+nA0ToWHoF9e9P4etEndFjv4Eh8xD4SNsf2
oIIOEWYbzPuxhbqJMfSMDiuBL2kz2NHX+POcq+nwNO2WLs9PPG/0SISuPqfD
rtDfIVRsZnV5u8SXON4P7Yw2YPv26tBevqbDBc58k53YmlsM1w6+pYPn8N0I
M+zZCNP+8Xd0WNyzaGCF7aBhlyFQTwcvFWLTeWy5KNeA/AY6bBjrsfXH7uv3
sdJtosP2FC+hMOzD0dfFz36kg5Dd5p/J2AKDyUt/WugwsjFOLAu7WetuR+wn
OqTuE36djx1+/X6FXBsdDnMPjBdj7xx6kvSsnQ5qvYIx5dir2jVe+7/S4YjV
x9Iq7KqYRouhDjpo79O0qMH2Hm5X9/lGh6QMF69a7M06/QLc3TgeZ5KFG7Gn
Yydm7/TQ4Wpv95aP2IUjCy2afXSIdLL7/Anbbtvfh439dNgad3i2HVs6nj32
1CAdepe6Yzqwu0cFXf8bosP7C7T737CTdaX2RY7g+M4c2NuNfSBBceO6MTq0
+aed7cHmGd/M+XicDiLGG1h6sRu2b5/YM0GHsSIV2j+H3thd3zNJh9Mcna/+
/b/+d/OC8z/oMLtwrP/f/YgdVtdYf9LBUGAsqAv76c2zZ9Nn6OBEL73bie05
4W6oNofzm3Pa8Cu2ir6fXO08rh9Vn53asCcSQ5mP/0cHAcMU/lbsvMnYwZ+L
dDja7K3TjH0apdaE/KbDyfjckTpsanJOltgffP3gScG32B1TxYHFyzg/FQfe
vMC+Cc+sDRh0+K/zzK+n2Ga33mzvWKHDo6tKxWXYXNPN1HOrdEircR+/h/3e
oGN5bQ3HQ88rLwf7fziHvN4=
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  AxesLabel->{
    FormBox["a", TraditionalForm], 
    FormBox["J", TraditionalForm]},
  AxesOrigin->{2., 1.},
  PlotRange->{{1, 10}, {0.8287248875490171, 2.5316571472731213`}},
  PlotRangeClipping->True,
  PlotRangePadding->{
    Scaled[0.02], 
    Scaled[0.02]}]], "Output",
 CellChangeTimes->{3.5043511000401287`*^9}],

Cell[BoxData["\<\"conv_diff_J.pdf\"\>"], "Output",
 CellChangeTimes->{3.5043511004285994`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"J", "[", "a", "]"}]], "Input",
 CellChangeTimes->{{3.502246270259112*^9, 3.502246271873864*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"5", " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"a", " ", 
      RowBox[{"(", 
       RowBox[{"45", "+", 
        RowBox[{"2", " ", 
         SuperscriptBox["a", "4"]}]}], ")"}]}], "+", 
     RowBox[{"45", " ", "a", " ", 
      RowBox[{"Cos", "[", 
       RowBox[{"2", " ", "a"}], "]"}]}], "+", 
     RowBox[{"15", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", "3"}], "+", 
        SuperscriptBox["a", "2"]}], ")"}], " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"2", " ", "a"}], "]"}]}]}], ")"}]}], 
  RowBox[{"6", " ", 
   SuperscriptBox["a", "5"]}]]], "Output",
 CellChangeTimes->{3.50435110054865*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"pdf", "=", 
   RowBox[{"1", "/", 
    RowBox[{"(", 
     RowBox[{"amax", "-", "amin"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"mean", "=", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"J", "[", "a", "]"}], "*", "pdf"}], ",", 
    RowBox[{"{", 
     RowBox[{"a", ",", "amin", ",", "amax"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"N", "[", 
  RowBox[{"mean", ",", "20"}], "]"}]}], "Input",
 CellChangeTimes->{{3.5022406543285637`*^9, 3.502240669686335*^9}, {
   3.50224092462926*^9, 3.502240928636731*^9}, {3.503149771667551*^9, 
   3.50314977474372*^9}, {3.5031500954647617`*^9, 3.50315009608607*^9}, {
   3.504350966138465*^9, 3.5043509915388117`*^9}, {3.504351033365028*^9, 
   3.504351039972763*^9}, 3.504351089329297*^9}],

Cell[BoxData[
 FractionBox[
  RowBox[{"87960", "+", 
   RowBox[{"20000", " ", 
    RowBox[{"Cos", "[", "2", "]"}]}], "-", 
   RowBox[{"20", " ", 
    RowBox[{"Cos", "[", "20", "]"}]}], "-", 
   RowBox[{"30000", " ", 
    RowBox[{"Sin", "[", "2", "]"}]}], "+", 
   RowBox[{"3", " ", 
    RowBox[{"Sin", "[", "20", "]"}]}]}], "28800"]], "Output",
 CellChangeTimes->{3.5043511046049347`*^9}],

Cell[BoxData["1.8178026964862019597711778602557428844227167767573047991676`20.\
"], "Output",
 CellChangeTimes->{3.504351104606147*^9}]
}, Open  ]]
},
WindowSize->{640, 756},
WindowMargins->{{108, Automatic}, {Automatic, 0}},
FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (November 6, \
2010)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 3291, 80, 193, "Input"],
Cell[3861, 104, 490, 15, 27, "Output"],
Cell[4354, 121, 698, 23, 33, "Output"],
Cell[5055, 146, 50941, 851, 243, "Output"],
Cell[55999, 999, 92, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[56128, 1005, 702, 19, 27, "Input"],
Cell[56833, 1026, 1333, 42, 63, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[58203, 1073, 1745, 36, 58, "Input"],
Cell[59951, 1111, 15301, 258, 247, "Output"],
Cell[75255, 1371, 94, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[75386, 1377, 120, 2, 27, "Input"],
Cell[75509, 1381, 668, 22, 53, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[76214, 1408, 813, 20, 58, "Input"],
Cell[77030, 1430, 388, 11, 47, "Output"],
Cell[77421, 1443, 135, 2, 27, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
